{
 "cells": [
  {
   "cell_type": "markdown",
   "id": "1df9294b",
   "metadata": {},
   "source": [
    "# MSRE isothermal temperature coefficient\n",
    "\n",
    "In this example we will calculate the isothermal temperature reactivity feedback coefficient of the MSRE. To do that we will parametrize the fuel salt temperature and run 3 different cases (in reality you would want to run many more). We will then fit the data with a linear function and approximate the average coefficient to the slope of the curve. We will verify that the coefficient is negative and is approximately equal to $-10\\sim15 pcm/K$, as reported here:  [ORNL, pag.39](http://moltensalt.org/references/static/downloads/pdf/ORNL-4233.pdf) "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "id": "49285ef9",
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "rm: cannot remove '*.xml': No such file or directory\n",
      "rm: cannot remove '*.h5': No such file or directory\n"
     ]
    },
    {
     "data": {
      "text/plain": [
       "256"
      ]
     },
     "execution_count": 1,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "%matplotlib inline\n",
    "\n",
    "import openmc\n",
    "import numpy as np\n",
    "from math import log10\n",
    "import matplotlib.pyplot as plt\n",
    "import os \n",
    "\n",
    "os.system('rm *.xml *.h5')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "id": "6d36f317",
   "metadata": {},
   "outputs": [],
   "source": [
    "# Set the same temperature to fuel salt, moderator and inor vessel\n",
    "def build_materials(temp):\n",
    "    # Fuel salt\n",
    "    salt_density = 2.3275 * (1 -1.18e-4 * ((temp*9/5+32) - 1181))\n",
    "    salt = openmc.Material(name=\"salt\", temperature = temp)\n",
    "    salt.add_nuclide('Li6',1.31480070E-05)\n",
    "    salt.add_nuclide('Li7', 0.262960140146177)\n",
    "    salt.add_nuclide('Be9',1.1863E-01)\n",
    "    salt.add_nuclide('Zr90',1.0543E-02)\n",
    "    salt.add_nuclide('Zr91',2.2991E-03)\n",
    "    salt.add_nuclide('Zr92',3.5142E-03)\n",
    "    salt.add_nuclide('Zr94',3.5613E-03)\n",
    "    salt.add_nuclide('Zr96',5.7375E-04)\n",
    "    salt.add_nuclide('Hf174',8.3786E-10)\n",
    "    salt.add_nuclide('Hf176',2.7545E-08)\n",
    "    salt.add_nuclide('Hf177',9.7401E-08)\n",
    "    salt.add_nuclide('Hf178',1.4285E-07)\n",
    "    salt.add_nuclide('Hf179',7.1323E-08)\n",
    "    salt.add_nuclide('Hf180',1.8370E-07)\n",
    "    salt.add_nuclide('U234',1.034276246E-05)\n",
    "    salt.add_nuclide('U235',1.009695816E-03)\n",
    "    salt.add_nuclide('U236',4.227809892E-06)\n",
    "    salt.add_nuclide('U238',2.168267822E-03)\n",
    "    salt.add_nuclide('Fe54',2.8551E-06)\n",
    "    salt.add_nuclide('Fe56',4.4818E-05)\n",
    "    salt.add_nuclide('Fe57',1.0350E-06)\n",
    "    salt.add_nuclide('Fe58',1.3775E-07)\n",
    "    salt.add_nuclide('Cr50',2.1224E-06)\n",
    "    salt.add_nuclide('Cr52',4.0928E-05)\n",
    "    salt.add_nuclide('Cr53',4.6409E-06)\n",
    "    salt.add_nuclide('Cr54',1.1552E-06)\n",
    "    salt.add_nuclide('Ni58',5.8597E-06)\n",
    "    salt.add_nuclide('Ni60',2.2571E-06)\n",
    "    salt.add_nuclide('Ni61',9.8117E-08)\n",
    "    salt.add_nuclide('Ni62',3.1284E-07)\n",
    "    salt.add_nuclide('Ni64',7.9671E-08)\n",
    "    salt.add_nuclide('O16',5.1437E-04)\n",
    "    salt.add_nuclide('O17',1.8927E-07)\n",
    "    salt.add_nuclide('O18',9.6440E-07)\n",
    "    salt.add_nuclide('F19',5.9409E-01)\n",
    "    salt.set_density('g/cm3', salt_density)\n",
    "\n",
    "    #Moderator graphite block\n",
    "    graphite = openmc.Material(name='graphite', temperature = temp)\n",
    "    graphite.set_density('g/cm3',1.86)\n",
    "    graphite.add_nuclide('C12',1)\n",
    "    graphite.add_s_alpha_beta('c_Graphite')\n",
    "\n",
    "    #inor-8\n",
    "    inor = openmc.Material(name='inor-8', temperature = temp)\n",
    "    inor.set_density('g/cm3',8.7745)\n",
    "    inor.add_element('Ni',68.5,'wo')\n",
    "    inor.add_element('Mo',16.5,'wo')\n",
    "    inor.add_element('Cr',7,'wo')\n",
    "    inor.add_element('Fe',5,'wo')\n",
    "    inor.add_element('C',0.06,'wo')\n",
    "    inor.add_element('Al',0.25,'wo')\n",
    "    inor.add_element('Ti',0.25,'wo')\n",
    "    inor.add_element('S',0.02,'wo')\n",
    "    inor.add_element('Mn',1.0,'wo')\n",
    "    inor.add_element('Si',1.0,'wo')\n",
    "    inor.add_element('Cu',0.35,'wo')\n",
    "    inor.add_element('B',0.010,'wo')\n",
    "    inor.add_element('W',0.5,'wo')\n",
    "    inor.add_element('P',0.015,'wo')\n",
    "    inor.add_element('Co',0.2,'wo')\n",
    "\n",
    "    #helium\n",
    "    helium = openmc.Material(name='helium')\n",
    "    helium.add_element('He',1.0)\n",
    "    helium.set_density('g/cm3',1.03*(10**-4))\n",
    "\n",
    "    #Control rods inconel clad\n",
    "    trace = 0.01\n",
    "    inconel = openmc.Material(name='inconel', temperature = 65.6 + 273.15)\n",
    "    inconel.add_element('Ni',78.5,percent_type='wo')\n",
    "    inconel.add_element('Cr',14.0,percent_type='wo')\n",
    "    inconel.add_element('Fe',6.5,percent_type='wo')\n",
    "    inconel.add_element('Mn',0.25,percent_type='wo')\n",
    "    inconel.add_element('Si',0.25,percent_type='wo')\n",
    "    inconel.add_element('Cu',0.2,percent_type='wo')\n",
    "    inconel.add_element('Co',0.2,percent_type='wo')\n",
    "    inconel.add_element('Al',0.2,percent_type='wo')\n",
    "    inconel.add_element('Ti',0.2,percent_type='wo')\n",
    "    inconel.add_element('Ta',0.5,percent_type='wo')\n",
    "    inconel.add_element('W',0.5,percent_type='wo')\n",
    "    inconel.add_element('Zn',0.2,percent_type='wo')\n",
    "    inconel.add_element('Zr',0.1,percent_type='wo')\n",
    "    inconel.add_element('C',trace,percent_type='wo')\n",
    "    inconel.add_element('Mo',trace,percent_type='wo')\n",
    "    inconel.add_element('Ag',trace,percent_type='wo')\n",
    "    inconel.add_element('B',trace,percent_type='wo')\n",
    "    inconel.add_element('Ba',trace,percent_type='wo')\n",
    "    inconel.add_element('Be',trace,percent_type='wo')\n",
    "    inconel.add_element('Ca',trace,percent_type='wo')\n",
    "    inconel.add_element('Cd',trace,percent_type='wo')\n",
    "    inconel.add_element('V',trace,percent_type='wo')\n",
    "    inconel.add_element('Sn',trace,percent_type='wo')\n",
    "    inconel.add_element('Mg',trace,percent_type='wo')\n",
    "    inconel.set_density('g/cm3',8.5)\n",
    "\n",
    "    # SS316 control rod flexible hose\n",
    "    ss316 = openmc.Material(name='ss316', temperature = 65.6 + 273.15)\n",
    "    ss316.add_element('C',0.026,'wo')\n",
    "    ss316.add_element('Si',0.37,'wo')\n",
    "    ss316.add_element('Mn',0.16,'wo')\n",
    "    ss316.add_element('Cr',16.55,'wo')\n",
    "    ss316.add_element('Cu',0.16,'wo')\n",
    "    ss316.add_element('Ni',10,'wo')\n",
    "    ss316.add_element('P',0.029,'wo')\n",
    "    ss316.add_element('S',0.027,'wo')\n",
    "    ss316.add_element('Mo',2.02,'wo')\n",
    "    ss316.add_element('N',0.036,'wo')\n",
    "    ss316.add_element('Fe',70.622,'wo')\n",
    "    ss316.set_density('g/cm3',7.99)\n",
    "\n",
    "    #Control rods bushing posion material\n",
    "    Gd2O3 = openmc.Material()\n",
    "    Gd2O3.add_element('Gd',2)\n",
    "    Gd2O3.add_element('O',3)\n",
    "    Gd2O3.set_density('g/cm3',7.41)\n",
    "    Al2O3 = openmc.Material()\n",
    "    Al2O3.add_element('Al',2)\n",
    "    Al2O3.add_element('O',3)\n",
    "    Al2O3.set_density('g/cm3',3.95)\n",
    "    bush = openmc.Material.mix_materials([Gd2O3,Al2O3],[0.7,0.3],'wo')\n",
    "    bush.name='bush'\n",
    "    bush.temperature = 65.6 +273.15\n",
    "\n",
    "    #Concrete block\n",
    "    concrete = openmc.Material(name='concrete')\n",
    "    concrete.add_element('H',0.005,'wo')\n",
    "    concrete.add_element('O',0.496,'wo')\n",
    "    concrete.add_element('Si',0.314,'wo')\n",
    "    concrete.add_element('Ca',0.083,'wo')\n",
    "    concrete.add_element('Na',0.017,'wo')\n",
    "    concrete.add_element('Mn',0.002,'wo')\n",
    "    concrete.add_element('Al',0.046,'wo')\n",
    "    concrete.add_element('S',0.001,'wo')\n",
    "    concrete.add_element('K',0.019,'wo')\n",
    "    concrete.add_element('Fe',0.012,'wo')\n",
    "    concrete.set_density('g/cm3',2.35)\n",
    "\n",
    "    #Thermal shielding as material mix of water and SS304 50-50 vo\n",
    "    water = openmc.Material()\n",
    "    water.add_element('H',2)\n",
    "    water.add_element('O',1)\n",
    "    water.set_density('g/cm3',0.997)\n",
    "\n",
    "    #stainless steel 304\n",
    "    ss304 =  openmc.Material()\n",
    "    ss304.add_element('C',0.08,'wo')\n",
    "    ss304.add_element('Mn',2,'wo')\n",
    "    ss304.add_element('P',0.045,'wo')\n",
    "    ss304.add_element('S',0.03,'wo')\n",
    "    ss304.add_element('Si',0.75,'wo')\n",
    "    ss304.add_element('Cr',19,'wo')\n",
    "    ss304.add_element('Ni',10,'wo')\n",
    "    ss304.add_element('N',0.1,'wo')\n",
    "    ss304.add_element('Fe',67.995, 'wo')\n",
    "    ss304.set_density('g/cm3',7.93)\n",
    "    shield = openmc.Material.mix_materials([water,ss304],[0.5,0.5],'vo')\n",
    "    shield.temperature = 32.2 + 273.15\n",
    "    shield.name='steelwater'\n",
    "\n",
    "    # \"Careytemp 1600\" by Philip Carey Manufacturing Compamy (Cincinnati) \n",
    "    insulation=openmc.Material(name='insulation')\n",
    "    insulation.add_element('Si',1)\n",
    "    insulation.add_element('O',2)\n",
    "    insulation.set_density('g/cm3',0.16) \n",
    "\n",
    "    # Sand water (not sure about this material)\n",
    "    sandwater=openmc.Material(name='sandwater')\n",
    "    sandwater.add_element('Fe',3)\n",
    "    sandwater.add_element('O',4)\n",
    "    sandwater.set_density('g/cm3',6)\n",
    "\n",
    "    #Vessel anular steel\n",
    "    steel = openmc.Material(name='steel')\n",
    "    steel.add_element('Fe',1)\n",
    "    steel.set_density('g/cm3',7.85)\n",
    "\n",
    "    return openmc.Materials([salt,graphite,inor,helium,inconel,shield,concrete,steel,ss316,sandwater,insulation,bush])\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "id": "5fa45be2",
   "metadata": {},
   "outputs": [],
   "source": [
    "# Import h5m files\n",
    "core_h5m = '../h5m/msre_reactor_1e-2.h5m'\n",
    "cr_h5m = '../h5m/msre_control_rod_1e-2.h5m'\n",
    "\n",
    "# Create DAGMC universes out of h5m files\n",
    "core = openmc.DAGMCUniverse(filename=core_h5m, auto_geom_ids=True, universe_id=1)\n",
    "cr = openmc.DAGMCUniverse(filename=cr_h5m, auto_geom_ids=True, universe_id=2)\n",
    "\n",
    "# Create regions\n",
    "core_region = core.bounding_region()\n",
    "cr1_region = cr.bounding_region(boundary_type='transmission', starting_id=20000)\n",
    "\n",
    "# Extend control rod region, to include upwards translations\n",
    "cr1_region = cr1_region | cr1_region.translate([0,0,150])\n",
    "\n",
    "# Create control rod region 2 and 3 as translated region of control rod 1\n",
    "offset = 10.163255\n",
    "cr2_region = cr1_region.translate([-offset,0,0])\n",
    "cr3_region = cr1_region.translate([-offset,offset,0])\n",
    "\n",
    "# Create openmc Cells \n",
    "core_cell = openmc.Cell(region=~(cr1_region | cr2_region | cr3_region) & core_region , fill=core)\n",
    "cr1_cell = openmc.Cell(name='CR1', region=cr1_region, fill=cr)\n",
    "cr2_cell = openmc.Cell(name='CR2', region=cr2_region, fill=cr)\n",
    "cr3_cell = openmc.Cell(name='CR3', region=cr3_region, fill=cr)\n",
    "    \n",
    "#translate control rods at top position (fully withdrawn)\n",
    "inch_to_cm = 2.54\n",
    "start_pos = 19.2 #inches\n",
    "top_pos = 51 #inches\n",
    "setattr(cr1_cell, 'translation', [0, 0, start_pos + top_pos*inch_to_cm])\n",
    "setattr(cr2_cell, 'translation', [-offset, 0, start_pos + top_pos*2.54])\n",
    "setattr(cr3_cell, 'translation', [-offset, offset, start_pos + top_pos*2.54])\n",
    "\n",
    "# Create openmc Geometry object\n",
    "geometry = openmc.Geometry([core_cell,cr1_cell,cr2_cell,cr3_cell])"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "id": "7bb7b060",
   "metadata": {},
   "outputs": [],
   "source": [
    "settings = openmc.Settings()\n",
    "settings.temperature = {'method':'interpolation','range':(293.15,923.15)}\n",
    "settings.batches = 50\n",
    "settings.inactive = 20\n",
    "settings.particles = 30000\n",
    "settings.photon_transport = False\n",
    "source_area = openmc.stats.Box([-100., -100., 0.],[ 100.,  100.,  200.],only_fissionable = True)\n",
    "settings.source = openmc.Source(space=source_area)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "id": "3b6184d5",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "                                %%%%%%%%%%%%%%%\n",
      "                           %%%%%%%%%%%%%%%%%%%%%%%%\n",
      "                        %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%\n",
      "                      %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%\n",
      "                    %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%\n",
      "                   %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%\n",
      "                                    %%%%%%%%%%%%%%%%%%%%%%%%\n",
      "                                     %%%%%%%%%%%%%%%%%%%%%%%%\n",
      "                 ###############      %%%%%%%%%%%%%%%%%%%%%%%%\n",
      "                ##################     %%%%%%%%%%%%%%%%%%%%%%%\n",
      "                ###################     %%%%%%%%%%%%%%%%%%%%%%%\n",
      "                ####################     %%%%%%%%%%%%%%%%%%%%%%\n",
      "                #####################     %%%%%%%%%%%%%%%%%%%%%\n",
      "                ######################     %%%%%%%%%%%%%%%%%%%%\n",
      "                #######################     %%%%%%%%%%%%%%%%%%\n",
      "                 #######################     %%%%%%%%%%%%%%%%%\n",
      "                 ######################     %%%%%%%%%%%%%%%%%\n",
      "                  ####################     %%%%%%%%%%%%%%%%%\n",
      "                    #################     %%%%%%%%%%%%%%%%%\n",
      "                     ###############     %%%%%%%%%%%%%%%%\n",
      "                       ############     %%%%%%%%%%%%%%%\n",
      "                          ########     %%%%%%%%%%%%%%\n",
      "                                      %%%%%%%%%%%\n",
      "\n",
      "                 | The OpenMC Monte Carlo Code\n",
      "       Copyright | 2011-2022 MIT, UChicago Argonne LLC, and contributors\n",
      "         License | https://docs.openmc.org/en/latest/license.html\n",
      "         Version | 0.13.2\n",
      "        Git SHA1 | 030f73a8690ed19e91806e46c8caf338d252e74a\n",
      "       Date/Time | 2023-01-17 13:43:17\n",
      "   MPI Processes | 1\n",
      "  OpenMP Threads | 56\n",
      "\n",
      " Reading settings XML file...\n",
      " Reading cross sections XML file...\n",
      " Reading materials XML file...\n",
      " Reading geometry XML file...\n",
      "Using the DOUBLE-DOWN interface to Embree.\n",
      "Loading file h5m/msre_reactor_1e-2.h5m\n",
      "Initializing the GeomQueryTool...\n",
      "Using faceting tolerance: 0.01\n",
      "Building acceleration data structures...\n",
      "Implicit Complement assumed to be Vacuum\n",
      "Using the DOUBLE-DOWN interface to Embree.\n",
      "Loading file h5m/msre_control_rod_1e-2.h5m\n",
      "Initializing the GeomQueryTool...\n",
      "Using faceting tolerance: 0.01\n",
      "Building acceleration data structures...\n",
      "Implicit Complement assumed to be Vacuum\n",
      " Reading Li6 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Li6.h5\n",
      " Reading Li7 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Li7.h5\n",
      " Reading Be9 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Be9.h5\n",
      " Reading Zr90 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Zr90.h5\n",
      " Reading Zr91 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Zr91.h5\n",
      " Reading Zr92 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Zr92.h5\n",
      " Reading Zr94 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Zr94.h5\n",
      " Reading Zr96 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Zr96.h5\n",
      " WARNING: Negative value(s) found on probability table for nuclide Zr96 at 250K\n",
      " WARNING: Negative value(s) found on probability table for nuclide Zr96 at 294K\n",
      " WARNING: Negative value(s) found on probability table for nuclide Zr96 at 600K\n",
      " WARNING: Negative value(s) found on probability table for nuclide Zr96 at 900K\n",
      " WARNING: Negative value(s) found on probability table for nuclide Zr96 at 1200K\n",
      " Reading Hf174 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Hf174.h5\n",
      " Reading Hf176 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Hf176.h5\n",
      " Reading Hf177 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Hf177.h5\n",
      " Reading Hf178 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Hf178.h5\n",
      " Reading Hf179 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Hf179.h5\n",
      " Reading Hf180 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Hf180.h5\n",
      " Reading U234 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/U234.h5\n",
      " Reading U235 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/U235.h5\n",
      " Reading U236 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/U236.h5\n",
      " Reading U238 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/U238.h5\n",
      " Reading Fe54 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Fe54.h5\n",
      " Reading Fe56 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Fe56.h5\n",
      " Reading Fe57 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Fe57.h5\n",
      " Reading Fe58 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Fe58.h5\n",
      " Reading Cr50 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Cr50.h5\n",
      " Reading Cr52 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Cr52.h5\n",
      " Reading Cr53 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Cr53.h5\n",
      " Reading Cr54 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Cr54.h5\n",
      " Reading Ni58 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Ni58.h5\n",
      " Reading Ni60 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Ni60.h5\n",
      " Reading Ni61 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Ni61.h5\n",
      " Reading Ni62 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Ni62.h5\n",
      " Reading Ni64 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Ni64.h5\n",
      " Reading O16 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/O16.h5\n",
      " Reading O17 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/O17.h5\n",
      " Reading O18 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/O18.h5\n",
      " Reading F19 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/F19.h5\n",
      " Reading C12 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/C12.h5\n",
      " Reading Mo100 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Mo100.h5\n",
      " Reading Mo92 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Mo92.h5\n",
      " Reading Mo94 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Mo94.h5\n",
      " Reading Mo95 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Mo95.h5\n",
      " Reading Mo96 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Mo96.h5\n",
      " Reading Mo97 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Mo97.h5\n",
      " Reading Mo98 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Mo98.h5\n",
      " Reading C13 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/C13.h5\n",
      " Reading Al27 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Al27.h5\n",
      " Reading Ti46 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Ti46.h5\n",
      " Reading Ti47 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Ti47.h5\n",
      " Reading Ti48 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Ti48.h5\n",
      " Reading Ti49 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Ti49.h5\n",
      " Reading Ti50 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Ti50.h5\n",
      " Reading S32 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/S32.h5\n",
      " Reading S33 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/S33.h5\n",
      " Reading S34 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/S34.h5\n",
      " Reading S36 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/S36.h5\n",
      " Reading Mn55 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Mn55.h5\n",
      " Reading Si28 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Si28.h5\n",
      " Reading Si29 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Si29.h5\n",
      " Reading Si30 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Si30.h5\n",
      " Reading Cu63 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Cu63.h5\n",
      " Reading Cu65 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Cu65.h5\n",
      " Reading B10 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/B10.h5\n",
      " Reading B11 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/B11.h5\n",
      " Reading W180 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/W180.h5\n",
      " Reading W182 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/W182.h5\n",
      " Reading W183 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/W183.h5\n",
      " Reading W184 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/W184.h5\n",
      " Reading W186 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/W186.h5\n",
      " Reading P31 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/P31.h5\n",
      " Reading Co59 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Co59.h5\n",
      " Reading He3 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/He3.h5\n",
      " Reading He4 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/He4.h5\n",
      " Reading Ta180 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Ta180.h5\n",
      " Reading Ta181 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Ta181.h5\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      " Reading Zn64 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Zn64.h5\n",
      " Reading Zn66 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Zn66.h5\n",
      " Reading Zn67 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Zn67.h5\n",
      " Reading Zn68 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Zn68.h5\n",
      " Reading Zn70 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Zn70.h5\n",
      " Reading Ag107 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Ag107.h5\n",
      " Reading Ag109 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Ag109.h5\n",
      " Reading Ba130 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Ba130.h5\n",
      " Reading Ba132 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Ba132.h5\n",
      " Reading Ba134 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Ba134.h5\n",
      " Reading Ba135 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Ba135.h5\n",
      " Reading Ba136 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Ba136.h5\n",
      " Reading Ba137 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Ba137.h5\n",
      " Reading Ba138 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Ba138.h5\n",
      " Reading Ca40 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Ca40.h5\n",
      " Reading Ca42 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Ca42.h5\n",
      " Reading Ca43 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Ca43.h5\n",
      " Reading Ca44 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Ca44.h5\n",
      " Reading Ca46 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Ca46.h5\n",
      " Reading Ca48 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Ca48.h5\n",
      " Reading Cd106 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Cd106.h5\n",
      " WARNING: Negative value(s) found on probability table for nuclide Cd106 at 250K\n",
      " WARNING: Negative value(s) found on probability table for nuclide Cd106 at 294K\n",
      " WARNING: Negative value(s) found on probability table for nuclide Cd106 at 600K\n",
      " WARNING: Negative value(s) found on probability table for nuclide Cd106 at 900K\n",
      " WARNING: Negative value(s) found on probability table for nuclide Cd106 at\n",
      "          1200K\n",
      " Reading Cd108 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Cd108.h5\n",
      " Reading Cd110 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Cd110.h5\n",
      " Reading Cd111 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Cd111.h5\n",
      " Reading Cd112 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Cd112.h5\n",
      " Reading Cd113 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Cd113.h5\n",
      " Reading Cd114 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Cd114.h5\n",
      " Reading Cd116 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Cd116.h5\n",
      " Reading V50 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/V50.h5\n",
      " Reading V51 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/V51.h5\n",
      " Reading Sn112 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Sn112.h5\n",
      " Reading Sn114 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Sn114.h5\n",
      " Reading Sn115 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Sn115.h5\n",
      " Reading Sn116 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Sn116.h5\n",
      " Reading Sn117 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Sn117.h5\n",
      " Reading Sn118 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Sn118.h5\n",
      " Reading Sn119 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Sn119.h5\n",
      " Reading Sn120 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Sn120.h5\n",
      " Reading Sn122 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Sn122.h5\n",
      " Reading Sn124 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Sn124.h5\n",
      " Reading Mg24 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Mg24.h5\n",
      " Reading Mg25 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Mg25.h5\n",
      " Reading Mg26 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Mg26.h5\n",
      " Reading N14 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/N14.h5\n",
      " Reading N15 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/N15.h5\n",
      " Reading Gd152 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Gd152.h5\n",
      " Reading Gd154 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Gd154.h5\n",
      " Reading Gd155 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Gd155.h5\n",
      " Reading Gd156 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Gd156.h5\n",
      " Reading Gd157 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Gd157.h5\n",
      " Reading Gd158 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Gd158.h5\n",
      " Reading Gd160 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Gd160.h5\n",
      " Reading H1 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/H1.h5\n",
      " Reading H2 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/H2.h5\n",
      " Reading Na23 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Na23.h5\n",
      " Reading K39 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/K39.h5\n",
      " Reading K40 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/K40.h5\n",
      " Reading K41 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/K41.h5\n",
      " Reading c_Graphite from\n",
      " /home/lorenzo/Documents/cross_sections/endfb80_hdf5/c_Graphite.h5\n",
      " Minimum neutron data temperature: 250 K\n",
      " Maximum neutron data temperature: 1200 K\n",
      " Preparing distributed cell instances...\n",
      " Reading plot XML file...\n",
      " Writing summary.h5 file...\n",
      " Maximum neutron transport energy: 20000000 eV for Li6\n",
      " Initializing source particles...\n",
      "\n",
      " ====================>     K EIGENVALUE SIMULATION     <====================\n",
      "\n",
      "  Bat./Gen.      k            Average k\n",
      "  =========   ========   ====================\n",
      "        1/1    0.59413\n",
      "        2/1    0.87144\n",
      "        3/1    0.95932\n",
      "        4/1    0.98230\n",
      "        5/1    0.98710\n",
      "        6/1    1.00177\n",
      "        7/1    1.00787\n",
      "        8/1    0.99209\n",
      "        9/1    0.99546\n",
      "       10/1    0.98464\n",
      "       11/1    1.01001\n",
      "       12/1    1.01107\n",
      " WARNING: No intersection found with DAGMC cell 48, material 3\n",
      "       13/1    1.00053\n",
      "       14/1    0.99278\n",
      "       15/1    0.99531\n",
      "       16/1    1.00262\n",
      "       17/1    0.98915\n",
      "       18/1    0.99553\n",
      "       19/1    0.98836\n",
      "       20/1    0.99625\n",
      "       21/1    0.99275\n",
      "       22/1    0.99348    0.99312 +/- 0.00037\n",
      "       23/1    0.98095    0.98906 +/- 0.00406\n",
      "       24/1    1.00824    0.99386 +/- 0.00559\n",
      "       25/1    1.00497    0.99608 +/- 0.00487\n",
      "       26/1    0.99537    0.99596 +/- 0.00397\n",
      "       27/1    0.99525    0.99586 +/- 0.00336\n",
      "       28/1    0.99927    0.99629 +/- 0.00294\n",
      "       29/1    0.99275    0.99589 +/- 0.00262\n",
      "       30/1    0.99491    0.99579 +/- 0.00235\n",
      "       31/1    1.00042    0.99622 +/- 0.00217\n",
      "       32/1    0.99247    0.99590 +/- 0.00200\n",
      "       33/1    0.99000    0.99545 +/- 0.00190\n",
      "       34/1    1.00013    0.99578 +/- 0.00179\n",
      "       35/1    1.00463    0.99637 +/- 0.00177\n",
      "       36/1    0.99759    0.99645 +/- 0.00165\n",
      "       37/1    1.00124    0.99673 +/- 0.00158\n",
      "       38/1    1.00066    0.99695 +/- 0.00150\n",
      "       39/1    0.98342    0.99624 +/- 0.00159\n",
      "       40/1    0.98485    0.99567 +/- 0.00161\n",
      "       41/1    1.00215    0.99598 +/- 0.00156\n",
      "       42/1    0.99562    0.99596 +/- 0.00149\n",
      "       43/1    1.00358    0.99629 +/- 0.00146\n",
      "       44/1    0.99845    0.99638 +/- 0.00140\n",
      "       45/1    0.98804    0.99605 +/- 0.00139\n",
      "       46/1    0.99894    0.99616 +/- 0.00134\n",
      "       47/1    0.99234    0.99602 +/- 0.00130\n",
      "       48/1    0.98830    0.99574 +/- 0.00128\n",
      "       49/1    1.00189    0.99595 +/- 0.00125\n",
      "       50/1    0.99920    0.99606 +/- 0.00121\n",
      " Creating state point statepoint.50.h5...\n",
      "\n",
      " =======================>     TIMING STATISTICS     <=======================\n",
      "\n",
      " Total time for initialization     = 1.1035e+02 seconds\n",
      "   Reading cross sections          = 1.5148e+01 seconds\n",
      " Total time in simulation          = 9.9007e+01 seconds\n",
      "   Time in transport only          = 9.8742e+01 seconds\n",
      "   Time in inactive batches        = 3.9581e+01 seconds\n",
      "   Time in active batches          = 5.9426e+01 seconds\n",
      "   Time synchronizing fission bank = 1.3993e-01 seconds\n",
      "     Sampling source sites         = 1.2650e-01 seconds\n",
      "     SEND/RECV source sites        = 1.3055e-02 seconds\n",
      "   Time accumulating tallies       = 2.7440e-05 seconds\n",
      "   Time writing statepoints        = 1.3125e-02 seconds\n",
      " Total time for finalization       = 3.5460e-06 seconds\n",
      " Total time elapsed                = 2.0983e+02 seconds\n",
      " Calculation Rate (inactive)       = 15158.7 particles/second\n",
      " Calculation Rate (active)         = 15144.8 particles/second\n",
      "\n",
      " ============================>     RESULTS     <============================\n",
      "\n",
      " k-effective (Collision)     = 0.99553 +/- 0.00123\n",
      " k-effective (Track-length)  = 0.99606 +/- 0.00121\n",
      " k-effective (Absorption)    = 0.99435 +/- 0.00114\n",
      " Combined k-effective        = 0.99510 +/- 0.00102\n",
      " Leakage Fraction            = 0.00001 +/- 0.00000\n",
      "\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "                                %%%%%%%%%%%%%%%\n",
      "                           %%%%%%%%%%%%%%%%%%%%%%%%\n",
      "                        %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%\n",
      "                      %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%\n",
      "                    %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%\n",
      "                   %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%\n",
      "                                    %%%%%%%%%%%%%%%%%%%%%%%%\n",
      "                                     %%%%%%%%%%%%%%%%%%%%%%%%\n",
      "                 ###############      %%%%%%%%%%%%%%%%%%%%%%%%\n",
      "                ##################     %%%%%%%%%%%%%%%%%%%%%%%\n",
      "                ###################     %%%%%%%%%%%%%%%%%%%%%%%\n",
      "                ####################     %%%%%%%%%%%%%%%%%%%%%%\n",
      "                #####################     %%%%%%%%%%%%%%%%%%%%%\n",
      "                ######################     %%%%%%%%%%%%%%%%%%%%\n",
      "                #######################     %%%%%%%%%%%%%%%%%%\n",
      "                 #######################     %%%%%%%%%%%%%%%%%\n",
      "                 ######################     %%%%%%%%%%%%%%%%%\n",
      "                  ####################     %%%%%%%%%%%%%%%%%\n",
      "                    #################     %%%%%%%%%%%%%%%%%\n",
      "                     ###############     %%%%%%%%%%%%%%%%\n",
      "                       ############     %%%%%%%%%%%%%%%\n",
      "                          ########     %%%%%%%%%%%%%%\n",
      "                                      %%%%%%%%%%%\n",
      "\n",
      "                 | The OpenMC Monte Carlo Code\n",
      "       Copyright | 2011-2022 MIT, UChicago Argonne LLC, and contributors\n",
      "         License | https://docs.openmc.org/en/latest/license.html\n",
      "         Version | 0.13.2\n",
      "        Git SHA1 | 030f73a8690ed19e91806e46c8caf338d252e74a\n",
      "       Date/Time | 2023-01-17 13:46:56\n",
      "   MPI Processes | 1\n",
      "  OpenMP Threads | 56\n",
      "\n",
      " Reading settings XML file...\n",
      " Reading cross sections XML file...\n",
      " Reading materials XML file...\n",
      " Reading geometry XML file...\n",
      "Using the DOUBLE-DOWN interface to Embree.\n",
      "Loading file h5m/msre_reactor_1e-2.h5m\n",
      "Initializing the GeomQueryTool...\n",
      "Using faceting tolerance: 0.01\n",
      "Building acceleration data structures...\n",
      "Implicit Complement assumed to be Vacuum\n",
      "Using the DOUBLE-DOWN interface to Embree.\n",
      "Loading file h5m/msre_control_rod_1e-2.h5m\n",
      "Initializing the GeomQueryTool...\n",
      "Using faceting tolerance: 0.01\n",
      "Building acceleration data structures...\n",
      "Implicit Complement assumed to be Vacuum\n",
      " Reading Li6 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Li6.h5\n",
      " Reading Li7 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Li7.h5\n",
      " Reading Be9 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Be9.h5\n",
      " Reading Zr90 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Zr90.h5\n",
      " Reading Zr91 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Zr91.h5\n",
      " Reading Zr92 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Zr92.h5\n",
      " Reading Zr94 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Zr94.h5\n",
      " Reading Zr96 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Zr96.h5\n",
      " WARNING: Negative value(s) found on probability table for nuclide Zr96 at 250K\n",
      " WARNING: Negative value(s) found on probability table for nuclide Zr96 at 294K\n",
      " WARNING: Negative value(s) found on probability table for nuclide Zr96 at 600K\n",
      " WARNING: Negative value(s) found on probability table for nuclide Zr96 at 900K\n",
      " WARNING: Negative value(s) found on probability table for nuclide Zr96 at 1200K\n",
      " Reading Hf174 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Hf174.h5\n",
      " Reading Hf176 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Hf176.h5\n",
      " Reading Hf177 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Hf177.h5\n",
      " Reading Hf178 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Hf178.h5\n",
      " Reading Hf179 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Hf179.h5\n",
      " Reading Hf180 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Hf180.h5\n",
      " Reading U234 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/U234.h5\n",
      " Reading U235 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/U235.h5\n",
      " Reading U236 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/U236.h5\n",
      " Reading U238 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/U238.h5\n",
      " Reading Fe54 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Fe54.h5\n",
      " Reading Fe56 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Fe56.h5\n",
      " Reading Fe57 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Fe57.h5\n",
      " Reading Fe58 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Fe58.h5\n",
      " Reading Cr50 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Cr50.h5\n",
      " Reading Cr52 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Cr52.h5\n",
      " Reading Cr53 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Cr53.h5\n",
      " Reading Cr54 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Cr54.h5\n",
      " Reading Ni58 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Ni58.h5\n",
      " Reading Ni60 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Ni60.h5\n",
      " Reading Ni61 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Ni61.h5\n",
      " Reading Ni62 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Ni62.h5\n",
      " Reading Ni64 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Ni64.h5\n",
      " Reading O16 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/O16.h5\n",
      " Reading O17 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/O17.h5\n",
      " Reading O18 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/O18.h5\n",
      " Reading F19 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/F19.h5\n",
      " Reading C12 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/C12.h5\n",
      " Reading Mo100 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Mo100.h5\n",
      " Reading Mo92 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Mo92.h5\n",
      " Reading Mo94 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Mo94.h5\n",
      " Reading Mo95 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Mo95.h5\n",
      " Reading Mo96 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Mo96.h5\n",
      " Reading Mo97 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Mo97.h5\n",
      " Reading Mo98 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Mo98.h5\n",
      " Reading C13 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/C13.h5\n",
      " Reading Al27 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Al27.h5\n",
      " Reading Ti46 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Ti46.h5\n",
      " Reading Ti47 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Ti47.h5\n",
      " Reading Ti48 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Ti48.h5\n",
      " Reading Ti49 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Ti49.h5\n",
      " Reading Ti50 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Ti50.h5\n",
      " Reading S32 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/S32.h5\n",
      " Reading S33 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/S33.h5\n",
      " Reading S34 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/S34.h5\n",
      " Reading S36 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/S36.h5\n",
      " Reading Mn55 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Mn55.h5\n",
      " Reading Si28 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Si28.h5\n",
      " Reading Si29 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Si29.h5\n",
      " Reading Si30 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Si30.h5\n",
      " Reading Cu63 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Cu63.h5\n",
      " Reading Cu65 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Cu65.h5\n",
      " Reading B10 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/B10.h5\n",
      " Reading B11 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/B11.h5\n",
      " Reading W180 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/W180.h5\n",
      " Reading W182 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/W182.h5\n",
      " Reading W183 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/W183.h5\n",
      " Reading W184 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/W184.h5\n",
      " Reading W186 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/W186.h5\n",
      " Reading P31 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/P31.h5\n",
      " Reading Co59 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Co59.h5\n",
      " Reading He3 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/He3.h5\n",
      " Reading He4 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/He4.h5\n",
      " Reading Ta180 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Ta180.h5\n",
      " Reading Ta181 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Ta181.h5\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      " Reading Zn64 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Zn64.h5\n",
      " Reading Zn66 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Zn66.h5\n",
      " Reading Zn67 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Zn67.h5\n",
      " Reading Zn68 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Zn68.h5\n",
      " Reading Zn70 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Zn70.h5\n",
      " Reading Ag107 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Ag107.h5\n",
      " Reading Ag109 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Ag109.h5\n",
      " Reading Ba130 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Ba130.h5\n",
      " Reading Ba132 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Ba132.h5\n",
      " Reading Ba134 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Ba134.h5\n",
      " Reading Ba135 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Ba135.h5\n",
      " Reading Ba136 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Ba136.h5\n",
      " Reading Ba137 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Ba137.h5\n",
      " Reading Ba138 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Ba138.h5\n",
      " Reading Ca40 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Ca40.h5\n",
      " Reading Ca42 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Ca42.h5\n",
      " Reading Ca43 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Ca43.h5\n",
      " Reading Ca44 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Ca44.h5\n",
      " Reading Ca46 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Ca46.h5\n",
      " Reading Ca48 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Ca48.h5\n",
      " Reading Cd106 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Cd106.h5\n",
      " WARNING: Negative value(s) found on probability table for nuclide Cd106 at 250K\n",
      " WARNING: Negative value(s) found on probability table for nuclide Cd106 at 294K\n",
      " WARNING: Negative value(s) found on probability table for nuclide Cd106 at 600K\n",
      " WARNING: Negative value(s) found on probability table for nuclide Cd106 at 900K\n",
      " WARNING: Negative value(s) found on probability table for nuclide Cd106 at\n",
      "          1200K\n",
      " Reading Cd108 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Cd108.h5\n",
      " Reading Cd110 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Cd110.h5\n",
      " Reading Cd111 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Cd111.h5\n",
      " Reading Cd112 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Cd112.h5\n",
      " Reading Cd113 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Cd113.h5\n",
      " Reading Cd114 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Cd114.h5\n",
      " Reading Cd116 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Cd116.h5\n",
      " Reading V50 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/V50.h5\n",
      " Reading V51 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/V51.h5\n",
      " Reading Sn112 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Sn112.h5\n",
      " Reading Sn114 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Sn114.h5\n",
      " Reading Sn115 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Sn115.h5\n",
      " Reading Sn116 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Sn116.h5\n",
      " Reading Sn117 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Sn117.h5\n",
      " Reading Sn118 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Sn118.h5\n",
      " Reading Sn119 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Sn119.h5\n",
      " Reading Sn120 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Sn120.h5\n",
      " Reading Sn122 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Sn122.h5\n",
      " Reading Sn124 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Sn124.h5\n",
      " Reading Mg24 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Mg24.h5\n",
      " Reading Mg25 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Mg25.h5\n",
      " Reading Mg26 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Mg26.h5\n",
      " Reading N14 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/N14.h5\n",
      " Reading N15 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/N15.h5\n",
      " Reading Gd152 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Gd152.h5\n",
      " Reading Gd154 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Gd154.h5\n",
      " Reading Gd155 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Gd155.h5\n",
      " Reading Gd156 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Gd156.h5\n",
      " Reading Gd157 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Gd157.h5\n",
      " Reading Gd158 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Gd158.h5\n",
      " Reading Gd160 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Gd160.h5\n",
      " Reading H1 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/H1.h5\n",
      " Reading H2 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/H2.h5\n",
      " Reading Na23 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Na23.h5\n",
      " Reading K39 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/K39.h5\n",
      " Reading K40 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/K40.h5\n",
      " Reading K41 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/K41.h5\n",
      " Reading c_Graphite from\n",
      " /home/lorenzo/Documents/cross_sections/endfb80_hdf5/c_Graphite.h5\n",
      " Minimum neutron data temperature: 250 K\n",
      " Maximum neutron data temperature: 1200 K\n",
      " Preparing distributed cell instances...\n",
      " Reading plot XML file...\n",
      " Writing summary.h5 file...\n",
      " Maximum neutron transport energy: 20000000 eV for Li6\n",
      " Initializing source particles...\n",
      "\n",
      " ====================>     K EIGENVALUE SIMULATION     <====================\n",
      "\n",
      "  Bat./Gen.      k            Average k\n",
      "  =========   ========   ====================\n",
      "        1/1    0.59985\n",
      "        2/1    0.86992\n",
      "        3/1    0.93662\n",
      "        4/1    0.97317\n",
      "        5/1    0.98089\n",
      "        6/1    0.99006\n",
      "        7/1    0.97465\n",
      "        8/1    0.98408\n",
      "        9/1    0.97688\n",
      "       10/1    0.98133\n",
      "       11/1    0.99735\n",
      "       12/1    0.97989\n",
      "       13/1    0.98175\n",
      "       14/1    0.99775\n",
      "       15/1    0.99462\n",
      "       16/1    0.98558\n",
      "       17/1    0.98992\n",
      "       18/1    0.97169\n",
      "       19/1    0.98839\n",
      "       20/1    0.97831\n",
      "       21/1    0.99006\n",
      "       22/1    0.99916    0.99461 +/- 0.00455\n",
      "       23/1    0.99631    0.99518 +/- 0.00269\n",
      "       24/1    0.98196    0.99187 +/- 0.00381\n",
      "       25/1    0.99825    0.99315 +/- 0.00322\n",
      "       26/1    0.97554    0.99021 +/- 0.00394\n",
      "       27/1    0.97223    0.98764 +/- 0.00420\n",
      "       28/1    0.97858    0.98651 +/- 0.00381\n",
      "       29/1    0.98642    0.98650 +/- 0.00336\n",
      "       30/1    0.98457    0.98631 +/- 0.00301\n",
      "       31/1    0.97973    0.98571 +/- 0.00279\n",
      "       32/1    0.98574    0.98571 +/- 0.00255\n",
      "       33/1    0.98439    0.98561 +/- 0.00235\n",
      "       34/1    0.98536    0.98559 +/- 0.00217\n",
      "       35/1    0.98966    0.98587 +/- 0.00204\n",
      "       36/1    0.99909    0.98669 +/- 0.00208\n",
      "       37/1    0.99175    0.98699 +/- 0.00198\n",
      "       38/1    0.97021    0.98606 +/- 0.00208\n",
      "       39/1    0.98010    0.98574 +/- 0.00200\n",
      "       40/1    0.99120    0.98602 +/- 0.00191\n",
      "       41/1    0.97703    0.98559 +/- 0.00187\n",
      "       42/1    0.98874    0.98573 +/- 0.00179\n",
      "       43/1    1.00000    0.98635 +/- 0.00182\n",
      "       44/1    0.99398    0.98667 +/- 0.00177\n",
      "       45/1    0.99704    0.98709 +/- 0.00175\n",
      "       46/1    0.99588    0.98742 +/- 0.00171\n",
      "       47/1    0.99577    0.98773 +/- 0.00168\n",
      "       48/1    0.98891    0.98777 +/- 0.00162\n",
      "       49/1    0.99467    0.98801 +/- 0.00158\n",
      "       50/1    0.99553    0.98826 +/- 0.00154\n",
      " Creating state point statepoint.50.h5...\n",
      "\n",
      " =======================>     TIMING STATISTICS     <=======================\n",
      "\n",
      " Total time for initialization     = 1.0659e+02 seconds\n",
      "   Reading cross sections          = 1.5197e+01 seconds\n",
      " Total time in simulation          = 1.0105e+02 seconds\n",
      "   Time in transport only          = 1.0081e+02 seconds\n",
      "   Time in inactive batches        = 3.9431e+01 seconds\n",
      "   Time in active batches          = 6.1616e+01 seconds\n",
      "   Time synchronizing fission bank = 1.2863e-01 seconds\n",
      "     Sampling source sites         = 1.1769e-01 seconds\n",
      "     SEND/RECV source sites        = 1.0570e-02 seconds\n",
      "   Time accumulating tallies       = 2.3695e-05 seconds\n",
      "   Time writing statepoints        = 1.3315e-02 seconds\n",
      " Total time for finalization       = 3.4400e-06 seconds\n",
      " Total time elapsed                = 2.0811e+02 seconds\n",
      " Calculation Rate (inactive)       = 15216.6 particles/second\n",
      " Calculation Rate (active)         = 14606.6 particles/second\n",
      "\n",
      " ============================>     RESULTS     <============================\n",
      "\n",
      " k-effective (Collision)     = 0.98752 +/- 0.00155\n",
      " k-effective (Track-length)  = 0.98826 +/- 0.00154\n",
      " k-effective (Absorption)    = 0.98970 +/- 0.00114\n",
      " Combined k-effective        = 0.98972 +/- 0.00125\n",
      " Leakage Fraction            = 0.00001 +/- 0.00000\n",
      "\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "                                %%%%%%%%%%%%%%%\n",
      "                           %%%%%%%%%%%%%%%%%%%%%%%%\n",
      "                        %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%\n",
      "                      %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%\n",
      "                    %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%\n",
      "                   %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%\n",
      "                                    %%%%%%%%%%%%%%%%%%%%%%%%\n",
      "                                     %%%%%%%%%%%%%%%%%%%%%%%%\n",
      "                 ###############      %%%%%%%%%%%%%%%%%%%%%%%%\n",
      "                ##################     %%%%%%%%%%%%%%%%%%%%%%%\n",
      "                ###################     %%%%%%%%%%%%%%%%%%%%%%%\n",
      "                ####################     %%%%%%%%%%%%%%%%%%%%%%\n",
      "                #####################     %%%%%%%%%%%%%%%%%%%%%\n",
      "                ######################     %%%%%%%%%%%%%%%%%%%%\n",
      "                #######################     %%%%%%%%%%%%%%%%%%\n",
      "                 #######################     %%%%%%%%%%%%%%%%%\n",
      "                 ######################     %%%%%%%%%%%%%%%%%\n",
      "                  ####################     %%%%%%%%%%%%%%%%%\n",
      "                    #################     %%%%%%%%%%%%%%%%%\n",
      "                     ###############     %%%%%%%%%%%%%%%%\n",
      "                       ############     %%%%%%%%%%%%%%%\n",
      "                          ########     %%%%%%%%%%%%%%\n",
      "                                      %%%%%%%%%%%\n",
      "\n",
      "                 | The OpenMC Monte Carlo Code\n",
      "       Copyright | 2011-2022 MIT, UChicago Argonne LLC, and contributors\n",
      "         License | https://docs.openmc.org/en/latest/license.html\n",
      "         Version | 0.13.2\n",
      "        Git SHA1 | 030f73a8690ed19e91806e46c8caf338d252e74a\n",
      "       Date/Time | 2023-01-17 13:50:33\n",
      "   MPI Processes | 1\n",
      "  OpenMP Threads | 56\n",
      "\n",
      " Reading settings XML file...\n",
      " Reading cross sections XML file...\n",
      " Reading materials XML file...\n",
      " Reading geometry XML file...\n",
      "Using the DOUBLE-DOWN interface to Embree.\n",
      "Loading file h5m/msre_reactor_1e-2.h5m\n",
      "Initializing the GeomQueryTool...\n",
      "Using faceting tolerance: 0.01\n",
      "Building acceleration data structures...\n",
      "Implicit Complement assumed to be Vacuum\n",
      "Using the DOUBLE-DOWN interface to Embree.\n",
      "Loading file h5m/msre_control_rod_1e-2.h5m\n",
      "Initializing the GeomQueryTool...\n",
      "Using faceting tolerance: 0.01\n",
      "Building acceleration data structures...\n",
      "Implicit Complement assumed to be Vacuum\n",
      " Reading Li6 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Li6.h5\n",
      " Reading Li7 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Li7.h5\n",
      " Reading Be9 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Be9.h5\n",
      " Reading Zr90 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Zr90.h5\n",
      " Reading Zr91 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Zr91.h5\n",
      " Reading Zr92 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Zr92.h5\n",
      " Reading Zr94 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Zr94.h5\n",
      " Reading Zr96 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Zr96.h5\n",
      " WARNING: Negative value(s) found on probability table for nuclide Zr96 at 250K\n",
      " WARNING: Negative value(s) found on probability table for nuclide Zr96 at 294K\n",
      " WARNING: Negative value(s) found on probability table for nuclide Zr96 at 600K\n",
      " WARNING: Negative value(s) found on probability table for nuclide Zr96 at 900K\n",
      " WARNING: Negative value(s) found on probability table for nuclide Zr96 at 1200K\n",
      " Reading Hf174 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Hf174.h5\n",
      " Reading Hf176 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Hf176.h5\n",
      " Reading Hf177 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Hf177.h5\n",
      " Reading Hf178 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Hf178.h5\n",
      " Reading Hf179 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Hf179.h5\n",
      " Reading Hf180 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Hf180.h5\n",
      " Reading U234 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/U234.h5\n",
      " Reading U235 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/U235.h5\n",
      " Reading U236 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/U236.h5\n",
      " Reading U238 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/U238.h5\n",
      " Reading Fe54 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Fe54.h5\n",
      " Reading Fe56 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Fe56.h5\n",
      " Reading Fe57 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Fe57.h5\n",
      " Reading Fe58 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Fe58.h5\n",
      " Reading Cr50 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Cr50.h5\n",
      " Reading Cr52 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Cr52.h5\n",
      " Reading Cr53 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Cr53.h5\n",
      " Reading Cr54 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Cr54.h5\n",
      " Reading Ni58 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Ni58.h5\n",
      " Reading Ni60 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Ni60.h5\n",
      " Reading Ni61 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Ni61.h5\n",
      " Reading Ni62 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Ni62.h5\n",
      " Reading Ni64 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Ni64.h5\n",
      " Reading O16 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/O16.h5\n",
      " Reading O17 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/O17.h5\n",
      " Reading O18 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/O18.h5\n",
      " Reading F19 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/F19.h5\n",
      " Reading C12 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/C12.h5\n",
      " Reading Mo100 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Mo100.h5\n",
      " Reading Mo92 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Mo92.h5\n",
      " Reading Mo94 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Mo94.h5\n",
      " Reading Mo95 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Mo95.h5\n",
      " Reading Mo96 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Mo96.h5\n",
      " Reading Mo97 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Mo97.h5\n",
      " Reading Mo98 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Mo98.h5\n",
      " Reading C13 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/C13.h5\n",
      " Reading Al27 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Al27.h5\n",
      " Reading Ti46 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Ti46.h5\n",
      " Reading Ti47 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Ti47.h5\n",
      " Reading Ti48 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Ti48.h5\n",
      " Reading Ti49 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Ti49.h5\n",
      " Reading Ti50 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Ti50.h5\n",
      " Reading S32 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/S32.h5\n",
      " Reading S33 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/S33.h5\n",
      " Reading S34 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/S34.h5\n",
      " Reading S36 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/S36.h5\n",
      " Reading Mn55 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Mn55.h5\n",
      " Reading Si28 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Si28.h5\n",
      " Reading Si29 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Si29.h5\n",
      " Reading Si30 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Si30.h5\n",
      " Reading Cu63 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Cu63.h5\n",
      " Reading Cu65 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Cu65.h5\n",
      " Reading B10 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/B10.h5\n",
      " Reading B11 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/B11.h5\n",
      " Reading W180 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/W180.h5\n",
      " Reading W182 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/W182.h5\n",
      " Reading W183 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/W183.h5\n",
      " Reading W184 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/W184.h5\n",
      " Reading W186 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/W186.h5\n",
      " Reading P31 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/P31.h5\n",
      " Reading Co59 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Co59.h5\n",
      " Reading He3 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/He3.h5\n",
      " Reading He4 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/He4.h5\n",
      " Reading Ta180 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Ta180.h5\n",
      " Reading Ta181 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Ta181.h5\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      " Reading Zn64 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Zn64.h5\n",
      " Reading Zn66 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Zn66.h5\n",
      " Reading Zn67 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Zn67.h5\n",
      " Reading Zn68 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Zn68.h5\n",
      " Reading Zn70 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Zn70.h5\n",
      " Reading Ag107 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Ag107.h5\n",
      " Reading Ag109 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Ag109.h5\n",
      " Reading Ba130 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Ba130.h5\n",
      " Reading Ba132 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Ba132.h5\n",
      " Reading Ba134 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Ba134.h5\n",
      " Reading Ba135 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Ba135.h5\n",
      " Reading Ba136 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Ba136.h5\n",
      " Reading Ba137 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Ba137.h5\n",
      " Reading Ba138 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Ba138.h5\n",
      " Reading Ca40 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Ca40.h5\n",
      " Reading Ca42 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Ca42.h5\n",
      " Reading Ca43 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Ca43.h5\n",
      " Reading Ca44 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Ca44.h5\n",
      " Reading Ca46 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Ca46.h5\n",
      " Reading Ca48 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Ca48.h5\n",
      " Reading Cd106 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Cd106.h5\n",
      " WARNING: Negative value(s) found on probability table for nuclide Cd106 at 250K\n",
      " WARNING: Negative value(s) found on probability table for nuclide Cd106 at 294K\n",
      " WARNING: Negative value(s) found on probability table for nuclide Cd106 at 600K\n",
      " WARNING: Negative value(s) found on probability table for nuclide Cd106 at 900K\n",
      " WARNING: Negative value(s) found on probability table for nuclide Cd106 at\n",
      "          1200K\n",
      " Reading Cd108 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Cd108.h5\n",
      " Reading Cd110 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Cd110.h5\n",
      " Reading Cd111 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Cd111.h5\n",
      " Reading Cd112 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Cd112.h5\n",
      " Reading Cd113 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Cd113.h5\n",
      " Reading Cd114 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Cd114.h5\n",
      " Reading Cd116 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Cd116.h5\n",
      " Reading V50 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/V50.h5\n",
      " Reading V51 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/V51.h5\n",
      " Reading Sn112 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Sn112.h5\n",
      " Reading Sn114 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Sn114.h5\n",
      " Reading Sn115 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Sn115.h5\n",
      " Reading Sn116 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Sn116.h5\n",
      " Reading Sn117 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Sn117.h5\n",
      " Reading Sn118 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Sn118.h5\n",
      " Reading Sn119 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Sn119.h5\n",
      " Reading Sn120 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Sn120.h5\n",
      " Reading Sn122 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Sn122.h5\n",
      " Reading Sn124 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Sn124.h5\n",
      " Reading Mg24 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Mg24.h5\n",
      " Reading Mg25 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Mg25.h5\n",
      " Reading Mg26 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Mg26.h5\n",
      " Reading N14 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/N14.h5\n",
      " Reading N15 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/N15.h5\n",
      " Reading Gd152 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Gd152.h5\n",
      " Reading Gd154 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Gd154.h5\n",
      " Reading Gd155 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Gd155.h5\n",
      " Reading Gd156 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Gd156.h5\n",
      " Reading Gd157 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Gd157.h5\n",
      " Reading Gd158 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Gd158.h5\n",
      " Reading Gd160 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Gd160.h5\n",
      " Reading H1 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/H1.h5\n",
      " Reading H2 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/H2.h5\n",
      " Reading Na23 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/Na23.h5\n",
      " Reading K39 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/K39.h5\n",
      " Reading K40 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/K40.h5\n",
      " Reading K41 from /home/lorenzo/Documents/cross_sections/endfb80_hdf5/K41.h5\n",
      " Reading c_Graphite from\n",
      " /home/lorenzo/Documents/cross_sections/endfb80_hdf5/c_Graphite.h5\n",
      " Minimum neutron data temperature: 250 K\n",
      " Maximum neutron data temperature: 1200 K\n",
      " Preparing distributed cell instances...\n",
      " Reading plot XML file...\n",
      " Writing summary.h5 file...\n",
      " Maximum neutron transport energy: 20000000 eV for Li6\n",
      " Initializing source particles...\n",
      "\n",
      " ====================>     K EIGENVALUE SIMULATION     <====================\n",
      "\n",
      "  Bat./Gen.      k            Average k\n",
      "  =========   ========   ====================\n",
      "        1/1    0.58567\n",
      "        2/1    0.86522\n",
      "        3/1    0.95393\n",
      "        4/1    0.97631\n",
      "        5/1    0.96262\n",
      "        6/1    0.97669\n",
      "        7/1    0.96797\n",
      "        8/1    0.96948\n",
      "        9/1    0.97837\n",
      "       10/1    0.96869\n",
      "       11/1    0.96878\n",
      "       12/1    0.97342\n",
      "       13/1    0.98926\n",
      "       14/1    0.97921\n",
      "       15/1    0.97947\n",
      "       16/1    0.99110\n",
      "       17/1    0.99039\n",
      "       18/1    0.99421\n",
      "       19/1    0.97364\n",
      "       20/1    0.99937\n",
      "       21/1    0.99256\n",
      "       22/1    0.98223    0.98739 +/- 0.00516\n",
      "       23/1    0.98421    0.98633 +/- 0.00316\n",
      "       24/1    0.99580    0.98870 +/- 0.00326\n",
      "       25/1    0.97299    0.98556 +/- 0.00403\n",
      "       26/1    0.97458    0.98373 +/- 0.00376\n",
      "       27/1    0.97297    0.98219 +/- 0.00353\n",
      "       28/1    0.98161    0.98212 +/- 0.00306\n",
      "       29/1    0.97695    0.98154 +/- 0.00276\n",
      "       30/1    0.97210    0.98060 +/- 0.00264\n",
      "       31/1    0.99224    0.98166 +/- 0.00261\n",
      "       32/1    0.98297    0.98177 +/- 0.00239\n",
      "       33/1    0.97624    0.98134 +/- 0.00224\n",
      "       34/1    0.98156    0.98136 +/- 0.00207\n",
      "       35/1    0.98727    0.98175 +/- 0.00197\n",
      "       36/1    0.98541    0.98198 +/- 0.00186\n",
      "       37/1    0.98024    0.98188 +/- 0.00175\n",
      "       38/1    0.97897    0.98172 +/- 0.00165\n",
      "       39/1    0.98925    0.98211 +/- 0.00161\n",
      "       40/1    0.98232    0.98212 +/- 0.00153\n",
      "       41/1    0.97919    0.98198 +/- 0.00146\n",
      "       42/1    0.99294    0.98248 +/- 0.00148\n",
      "       43/1    0.98525    0.98260 +/- 0.00142\n",
      "       44/1    0.97696    0.98237 +/- 0.00138\n",
      "       45/1    0.99247    0.98277 +/- 0.00138\n",
      "       46/1    0.97144    0.98234 +/- 0.00140\n",
      "       47/1    0.97344    0.98201 +/- 0.00139\n",
      "       48/1    0.99184    0.98236 +/- 0.00138\n",
      "       49/1    0.98350    0.98240 +/- 0.00133\n",
      "       50/1    0.97664    0.98220 +/- 0.00130\n",
      " Creating state point statepoint.50.h5...\n",
      "\n",
      " =======================>     TIMING STATISTICS     <=======================\n",
      "\n",
      " Total time for initialization     = 1.0648e+02 seconds\n",
      "   Reading cross sections          = 1.4451e+01 seconds\n",
      " Total time in simulation          = 1.0340e+02 seconds\n",
      "   Time in transport only          = 1.0315e+02 seconds\n",
      "   Time in inactive batches        = 4.1904e+01 seconds\n",
      "   Time in active batches          = 6.1495e+01 seconds\n",
      "   Time synchronizing fission bank = 1.3368e-01 seconds\n",
      "     Sampling source sites         = 1.2047e-01 seconds\n",
      "     SEND/RECV source sites        = 1.2851e-02 seconds\n",
      "   Time accumulating tallies       = 2.2970e-05 seconds\n",
      "   Time writing statepoints        = 1.2931e-02 seconds\n",
      " Total time for finalization       = 4.8340e-06 seconds\n",
      " Total time elapsed                = 2.1040e+02 seconds\n",
      " Calculation Rate (inactive)       = 14318.3 particles/second\n",
      " Calculation Rate (active)         = 14635.4 particles/second\n",
      "\n",
      " ============================>     RESULTS     <============================\n",
      "\n",
      " k-effective (Collision)     = 0.98212 +/- 0.00133\n",
      " k-effective (Track-length)  = 0.98220 +/- 0.00130\n",
      " k-effective (Absorption)    = 0.98054 +/- 0.00135\n",
      " Combined k-effective        = 0.98143 +/- 0.00120\n",
      " Leakage Fraction            = 0.00001 +/- 0.00000\n",
      "\n"
     ]
    }
   ],
   "source": [
    "# Define temperature to simulate in Celsius\n",
    "temperatures = [598.9,648.9,698.9]\n",
    "keff = []\n",
    "\n",
    "for temp in temperatures:\n",
    "    model = openmc.model.Model(geometry, build_materials(temp+273.15), settings)\n",
    "    res = model.run()\n",
    "    with openmc.StatePoint(res) as sp:\n",
    "        keff.append((sp.keff.n, sp.keff.s))"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "f3b687f4",
   "metadata": {},
   "source": [
    "We will now fit the data with a linear function and approximate the slope of the curve to the average isothermal reactivity coefficient. To do that we will use the `scipy` library functionalities:\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "id": "658554f6",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "Text(0, 0.5, '$\\\\rho\\\\,\\\\pm\\\\sigma$')"
      ]
     },
     "execution_count": 6,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAlYAAAGwCAYAAABrUCsdAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjYuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8o6BhiAAAACXBIWXMAAA9hAAAPYQGoP6dpAACEK0lEQVR4nOzdd1gUV9vH8e8uHXFBFFhR7AWxCxYUOxFrjOJjibGXqKjRmER9kxhTjUmMvaYZo6aYGFuMBrErNpTYsYuKgIqASGfn/WPjmo00den357r2etjZc2bumfCwP8/MnFEpiqIghBBCCCGem7qgCxBCCCGEKC4kWAkhhBBCmIgEKyGEEEIIE5FgJYQQQghhIhKshBBCCCFMRIKVEEIIIYSJSLASQgghhDAR84IuoKTR6XRERERQunRpVCpVQZcjhBBCiFxQFIUHDx7g6uqKWp31uJQEq3wWERGBm5tbQZchhBBCiGdw48YNKlasmOXnEqzyWenSpQH9fxiNRlPA1QghhBAiN+Lj43FzczN8j2dFglU+e3T6T6PRSLASQgghipicLuORi9eFEEIIIUxEgpUQQgghhIlIsBJCCCGEMBG5xkoIIUShoNPpSE1NLegyRAllYWGBmZnZc69HgpUQQogCl5qaytWrV9HpdAVdiijBHBwc0Gq1zzXPpAQrIYQQBUpRFG7fvo2ZmRlubm7ZTr4oRF5QFIXExESio6MBKF++/DOvS4KVEEKIApWenk5iYiKurq7Y2toWdDmihLKxsQEgOjoaZ2fnZz4tKP8sEEIIUaAyMjIAsLS0LOBKREn3KNinpaU98zokWAkhhCgU5PmpoqCZ4ndQgpUQQgghhIlIsBJCCCGEMBEJVkIIIYTI1O7du1GpVKhUKl566aWCLueZ5ed+FLlgFRMTw8CBA9FoNDg4ODBixAgSEhKy7ZOcnExAQABly5bFzs4Of39/oqKijNqEh4fTrVs3bG1tcXZ25s033yQ9PT3T9R04cABzc3MaNWpkqt0SQghRDKxfv55OnTpRtmxZVCoVoaGhWbZVFIUuXbqgUqnYsGFDtutVFIUZM2ZQvnx5bGxs8PX15eLFi0ZtcvP9ePLkSVq3bo21tTVubm589tlnudqvsLAwVq5cmau2+WHYsGG88847hve7du2ia9eulC1bFltbWzw8PJgyZQq3bt0CoGXLlty+fZu+ffvmeW1FLlgNHDiQM2fOEBgYyJYtW9i7dy+jR4/Ots/kyZPZvHkz69atY8+ePURERNC7d2/D5xkZGXTr1o3U1FQOHjzI999/z8qVK5kxY8YT64qNjWXw4MF07NjR5PsmhBCiaHv48CE+Pj7Mnj07x7bz5s3L9cXSn332GQsWLGDZsmUcPnyYUqVK4efnR3JysqFNTt+P8fHxdOrUicqVKxMSEsLnn3/OzJkzWbFiRY7bd3Z2xsHBIVe15rWMjAy2bNnCiy++CMDy5cvx9fVFq9Xy22+/cfbsWZYtW0ZcXBxz5swB9HecarVaw5QKeUopQs6ePasAytGjRw3L/vzzT0WlUim3bt3KtE9sbKxiYWGhrFu3zrDs3LlzCqAEBwcriqIoW7duVdRqtRIZGWlos3TpUkWj0SgpKSlG6+vXr5/yzjvvKO+9957SsGHDp96HuLg4BVDi4uKeum9WHqakKZWnblEqT92iPExJM9l6hRAiPyQlJSlnz55VkpKSjD9IScj6lfo0bRNz19aErl69qgDKiRMnMv38xIkTSoUKFZTbt28rgPL7779nuS6dTqdotVrl888/NyyLjY1VrKyslB9//FFRlNx9Py5ZskQpU6aM0ffa1KlTldq1a2e57V27dimAcv/+faPlbdu2VQICApSAgABFo9EoZcuWVd555x1Fp9MZ2iQnJytvvfWWUrFiRcXS0lKpXr268vXXXxutd9u2bUqjRo0Ua2trpX379kpUVJSydetWxd3dXSldurQyYMAA5eHDh0bb3rt3r1K+fHlFp9MpN27cUCwtLZVJkyZlWv9/6x4yZIjSs2fPLPc3y99FJfff30VqgtDg4GAcHBzw8vIyLPP19UWtVnP48GF69er1RJ+QkBDS0tLw9fU1LHN3d6dSpUoEBwfTokULgoODqV+/Pi4uLoY2fn5+jB07ljNnztC4cWMAvvvuO65cucLq1av56KOPclVzSkoKKSkphvfx8fFPvd9CCFEifeKa9Wc1O8HAdY/ff14D0hIzb1vZB4b98fj9vPqQeO/JdjPjnq3Op5SYmMjLL7/M4sWL0Wq1Oba/evUqkZGRRt9j9vb2NG/enODgYPr375+r78fg4GDatGljNF+Yn58fs2fP5v79+5QpU+ap9uP7779nxIgRHDlyhGPHjjF69GgqVarEqFGjABg8eDDBwcEsWLCAhg0bcvXqVe7evWu0jpkzZ7Jo0SJsbW3p27cvffv2xcrKirVr15KQkECvXr1YuHAhU6dONfTZtGkTPXr0QKVSsW7dOlJTU3nrrbcyrbEgRtmKVLCKjIzE2dnZaJm5uTmOjo5ERkZm2cfS0vKJg+vi4mLoExkZaRSqHn3+6DOAixcvMm3aNPbt24e5ee4P26xZs3j//fdz3V4IIUTxNnnyZFq2bEnPnj1z1f7R91Bm31P//h7L6fsxMjKSqlWrPrGOR589bbByc3Nj7ty5qFQqateuzalTp5g7dy6jRo3iwoUL/PLLLwQGBhoCYbVq1Z5Yx0cffUSrVq0AGDFiBNOnT+fy5cuGtn369GHXrl1GwWrjxo3MnTsX0H83azSa53oEjakVimA1bdq0HM9Hnzt3Lp+qeVJGRgYvv/wy77//PrVq1XqqvtOnT+f11183vI+Pj8fNzc3UJQohRPHzfxFZf6b6z+NG3ryUTdv/XE486dSz1/SPNWvW8Oqrrxre//nnn7Ru3TrHfps2bWLnzp2cOHHiuWsoaC1atDC6Rszb25s5c+aQkZFBaGgoZmZmtG3bNtt1NGjQwPCzi4sLtra2RgHMxcWFI0eOGN6fO3eOiIgIw3XOiqIUuollC0WwmjJlCkOHDs22TbVq1dBqtYYHJD6Snp5OTExMlsOpWq2W1NRUYmNjjUatoqKiDH20Wq3Rf7hHnz/67MGDBxw7dowTJ04wfvx4AHQ6HYqiYG5uzl9//UWHDh0y3b6VlRVWVlbZ7psQQohMWJYq+LZZePHFF2nevLnhfYUKFXLVb+fOnVy+fPmJsyj+/v60bt2a3bt3P9Hn0XdVVFSU0chMVFSU4e703Hw/arXaJ+6I//d3nSnl9iJxCwsLw88qlcro/aNlOp3O8H7Tpk288MILWFtbA1CrVi3i4uK4fft2oRm1KhR3BTo5OeHu7p7ty9LSEm9vb2JjYwkJCTH03blzJzqdzugX/N88PT2xsLAgKCjIsCwsLIzw8HC8vb0Bfco+deqU0S9lYGAgGo0GDw8PNBoNp06dIjQ01PAaM2YMtWvXJjQ0NMttCyGEKJ5Kly5NjRo1DK/cBolp06Zx8uRJo+8TgLlz5/Ldd99l2qdq1apotVqj77H4+HgOHz5s9D2W0/ejt7c3e/fuNXoOXmBgILVr137q04AAhw8fNnp/6NAhatasiZmZGfXr10en07Fnz56nXm92Nm7caHQKtU+fPlhaWmY5bURsbKxJt58bhWLEKrfq1KlD586dGTVqFMuWLSMtLY3x48fTv39/XF31FzneunWLjh07smrVKpo1a4a9vT0jRozg9ddfx9HREY1Gw4QJE/D29qZFixYAdOrUCQ8PDwYNGsRnn31GZGQk77zzDgEBAYbRpnr16hnV4uzsjLW19RPLhRBClFwxMTGEh4cTEaE/jRkWFgboR4T+/fqvSpUqGV3/5O7uzqxZs+jVqxcqlYpJkybx0UcfUbNmTapWrcq7776Lq6urYbLL3Hw/PrqkZcSIEUydOpXTp08zf/58w/VKTys8PJzXX3+dV199lePHj7Nw4ULD9AZVqlRhyJAhDB8+3HDx+vXr14mOjn7muaSio6M5duwYmzZtMix7dJ3X+PHjiY+PZ/DgwVSpUoWbN2+yatUq7OzsDDXllyIVrEB/Xnv8+PF07NgRtVqNv78/CxYsMHyelpZGWFgYiYmP7w6ZO3euoW1KSgp+fn4sWbLE8LmZmRlbtmxh7NixeHt7U6pUKYYMGcIHH3yQr/smhBCiaNu0aRPDhg0zvO/fvz8A7733HjNnzsz1esLCwoiLe3yX4ltvvcXDhw8ZPXo0sbGx+Pj4sG3bNsMpMcj5+9He3p6//vqLgIAAPD09KVeuHDNmzMhxLsisDB48mKSkJJo1a4aZmRmvvfaa0bqWLl3K//3f/zFu3Dju3btHpUqV+L//+79n2hbA5s2badasGeXKlTNaPm7cOGrVqsUXX3xBr169SEpKokqVKnTv3t3oGuf8olIURcn3rZZg8fHx2NvbExcXh0ajMck6E1PT8ZixHYCzH/hha1nk8rIQogRLTk7m6tWrVK1a1SgoiIK3e/du2rdvz/37942uC2vXrh2NGjVi3rx5+VbLiy++iI+PT5ZTK+TG0KFDiY2NzXKm++x+F3P7/V0orrESplFDdRPSU3JuKIQQQjyFihUrMmDAgAKtwcfH55lr2LdvH3Z2dqxZs8bEVT1JhjaKidIkstpyFtbffgM9F4Fbs4IuSQghRBHXvHlzwzMJ7ezsCrSW5xmp8vLyMtwokNf7IcGqmKiquo0ZGajvhsE3naDZaOj4LliVLujShBBCFFE2NjbUqFHjieWZTQtRmGW1H3lBTgUWEyeV6ryQ8jnpDQYAChxZDku84WJgQZcmhBBClBgSrIqRWEqT2n0RDPodHCpB3A1Y0wd+GwkZaTmvQAghhBDPRYJVcVS9A4w7BN7jHz/Kwcwi+z5CCCGEeG5yjVVxZVkK/D6Gur2hTOXHyx9EQUaKfkRLCCGEECYlI1bFXUVPKPXPZGqKAn+8DotbwKFloMso2NqEEMKEElPTqTLtD6pM+4PE1PSCLkeUUBKsSpLUh5B4D9Iewrap8K0fRJ8r6KqEEKLYUqlUWU5GaUpVqlQxyWSdplrP01q5cuUTD6YuqiRYlSRWdjB0K3SbA5al4eZRWNYads2SiUWFEOIp3blzh7Fjx1KpUiWsrKzQarX4+flx4MABQ5vbt2/TpUuXAqwyc1kFmaNHjz7zI27yW0GFwJzINVYljVoNTUdCrS7wxxS48Cfs+RTOboD+a6Fs9YKuUAghigR/f39SU1P5/vvvqVatGlFRUQQFBXHv3j1Dm8weuFyYOTk5FXQJRZ6MWBUD/37c42/Hb5KWocu5k30FGPAj9PkWbMtBWiKULlp/AIQQoqDExsayb98+Zs+eTfv27alcuTLNmjVj+vTpvPjii4Z2/z4VeO3aNVQqFb/88gutW7fGxsaGpk2bcuHCBY4ePYqXlxd2dnZ06dKFO3fuGNbRrl07Jk2aZLT9l156iaFDh2ZZ35dffkn9+vUpVaoUbm5ujBs3joSEBEA/ueewYcOIi4tDpVKhUqkMD4j+7yhQeHg4PXv2xM7ODo1GQ9++fYmKijJ8PnPmTBo1asQPP/xAlSpVsLe3p3///jx48CDb47dy5UoqVaqEra0tvXr1MgqjAJcvX6Znz564uLhgZ2dH06ZN2bFjh9ExuX79OpMnTzbsA8C9e/cYMGAAFSpUwNbWlvr16/Pjjz9mW4upSbAqBm7FJhl+fnfDGTrO2cMvx26QnlPAUqmgnj+MPwr91ujvJATQ6eB6cB5WLIQQWVMUhcTU9Gd6PfKs/f/9D9Xs2NnZYWdnx4YNG0hJebpLKd577z3eeecdjh8/jrm5OS+//DJvvfUW8+fPZ9++fVy6dIkZM2Y81Tr/S61Ws2DBAs6cOcP333/Pzp07DY+EadmyJfPmzUOj0XD79m1u377NG2+88cQ6dDodPXv2JCYmhj179hAYGMiVK1fo16+fUbvLly+zYcMGtmzZwpYtW9izZw+ffvpplrUdPnyYESNGMH78eEJDQ2nfvj0fffSRUZuEhAS6du1KUFAQJ06coHPnzvTo0YPw8HAA1q9fT8WKFfnggw8M+wD6hyh7enryxx9/cPr0aUaPHs2gQYM4cuTIcx3PpyGnAouBimVsDT+XLWVJeEwib/16ksW7LjG+fQ16Na6AuVk2GdrWUf965Ng3sPUNqN8XOn8KpcrmYfVCCGEsKS0Djxnbn2sdXh8FPVO/sx/4YWuZ81ejubk5K1euZNSoUSxbtowmTZrQtm1b+vfvT4MGDbLt+8Ybb+Dn5wfAa6+9xoABAwgKCqJVq1YAjBgxgpUrVz5T/Y/8e4SrSpUqfPTRR4wZM4YlS5ZgaWmJvb09KpUq21OVQUFBnDp1iqtXr+Lm5gbAqlWrqFu3LkePHqVp06aAPoCtXLmS0qX1j1AbNGgQQUFBfPzxx5mud/78+XTu3NkQ9GrVqsXBgwfZtm2boU3Dhg1p2LCh4f2HH37I77//zqZNmxg/fjyOjo6YmZlRunRpo32oUKGCUUicMGEC27dv55dffqFZs/x5hq6MWBUz2ye35u2udShbypLr9xJ589eTdPxyD7+G3Mx5BOuRh3f0E4ue+gUWN4WT6/RTNQghhDDw9/cnIiKCTZs20blzZ3bv3k2TJk1yDEX/Dl4uLi4A1K9f32hZdHT0c9W2Y8cOOnbsSIUKFShdujSDBg3i3r17JCYm5nod586dw83NzRCqADw8PHBwcODcucd3lFepUsUQqgDKly+fbf3nzp2jefPmRsu8vb2N3ickJPDGG29Qp04dHBwcsLOz49y5c4YRq6xkZGTw4YcfUr9+fRwdHbGzs2P79u059jMlGbEqZmwtzRnVphoDW1Ri9aHrLN9zhev3Enlj3d8s2nmRCR1q0rORa/YjWO3/D2r6waYJEH0G1o/Uh6xuX4KDW9b9hBDCBGwszDj7gd9T90tMTTeMVB17p2OuRp4y2/bTsLa25oUXXuCFF17g3XffZeTIkbz33nvZXv9kYfH4SRiPrg367zKd7vE/hNVq9ROnKNPSsn5M2bVr1+jevTtjx47l448/xtHRkf379zNixAhSU1OxtbXNsu+z+HftmdX/LN544w0CAwP54osvqFGjBjY2NvTp04fU1NRs+33++efMnz+fefPmGa4xmzRpUo79TElGrIopW0tzRrepzr6p7ZnexR3HUpZcu5fIlHV/88Lcvaw/nsMIVkVPGL0b2r8DZpZw8S9Y0gJOr8+3fRBClEwqlQpbS/Nnej3yrP0fBZ1n5eHhwcOHD5/3EBhxcnIyXEME+lGZ06dPZ9k+JCQEnU7HnDlzaNGiBbVq1SIiIsKojaWlJRkZ2U8SXadOHW7cuMGNGzcMy86ePUtsbCweHh7PuDf69R4+fNho2aFDh4zeHzhwgKFDh9KrVy/q16+PVqvl2rVrOe7DgQMH6NmzJ6+88goNGzakWrVqXLhw4ZlrfRYSrIo5W0tzXm1bnX1vtWfaPwHr6t2HvP7L33Sau5ffT9wkQ5fFaT5zS2j7JozZD24t9HcOOlbN3x0QQohC6N69e3To0IHVq1dz8uRJrl69yrp16/jss8/o2bOnSbfVoUMH/vjjD/744w/Onz/P2LFjiY2NzbJ9jRo1SEtLY+HChVy5coUffviBZcuWGbWpUqUKCQkJBAUFcffu3UxPEfr6+lK/fn0GDhzI8ePHOXLkCIMHD6Zt27Z4eXk98/5MnDiRbdu28cUXX3Dx4kUWLVpkdH0VQM2aNVm/fj2hoaH8/fffvPzyy0+MglWpUoW9e/dy69Yt7t69a+gXGBjIwYMHOXfuHK+++qrRXYz5QYJVCVHKypwx/wSsqZ3dKWNrwZW7D5n889+8MHcPG07cyjpgOdWGYX/C8O3g2vjx8mv7IT3/hleFEKKwsLOzo3nz5sydO5c2bdpQr1493n33XUaNGsWiRYtMuq3hw4czZMgQQ6ipVq0a7du3z7J9w4YN+fLLL5k9ezb16tVjzZo1zJo1y6hNy5YtGTNmDP369cPJyYnPPvvsifWoVCo2btxImTJlaNOmDb6+vlSrVo2ff/75ufanRYsWfPXVV8yfP5+GDRvy119/8c477xi1+fLLLylTpgwtW7akR48e+Pn50aRJE6M2H3zwAdeuXaN69eqG+bfeeecdmjRpgp+fH+3atUOr1fLSSy89V71PS6Xk9t5SYRLx8fHY29sTFxeHRqMxyToTU9MNd9Dk9o6WhJR0VgVfY8XeK8Qm6s/VV3MqxWsda9K9gStm6hyGw6PP6WdtL1sDXlwIbk2fez+EECVTcnIyV69epWrVqlhbWz/zep7lb6EQ/5bd72Juv79lxKqEsrMyZ1y7Guyf2oE3/WrjYGvBlTsPee2nUDrN3cPG0GxGsAAe3AZre7hzDr55Af6cBikJ+bcDQgghRCEkwaqEs7MyJ6B9Dfa91Z43/Wpjb2PB5X8Clt+8vWz6OyLzgFW9g35i0YYDAAUOL4Ul3nBxx5NthRAiH9hamnPt025c+7SbjFaJAiPBSgBQ2tqCgPY12D+1PW90qoW9jQWXohOY+OMJOs/by+a/I9D9N2DZOkKvZfDKenCoBHHhsMYfNr9WMDshhBBCFDAJVsJIaWsLxneoyb6p7ZnyQi001uZcjE5gwo8n6Dx/L1tOZhKwanSEcYegRYB+YtGyNQqmeCGEEKKASbASmdJYWzChY032T+vA6/8ErAtRCYxfqw9Yf5y8bRywLEtB509g9B5oPvbx8tsnIfbGkxsQQoj/kHupREEzxe+gBCuRLY21BRM71mTf1A5M9q1F6X8CVsDa43SZv4+tp/4TsMo3ALN/rm1IS4Zfh+knFj28Qv9wZyGE+A8zM/1s5/k5O7YQmXk0n9d/Z5N/GnJ1n8gVexsLXvOtydBWVfjuwFW+2X+VsKgHjFtzHHdtaV7rWBO/ulrU/56mIek+2JaDe5fgzzfh1Dr91AzO7gW3I0KIQsfc3BxbW1vu3LmDhYUFarX8m1/kL0VRSExMJDo6GgcHB0PYfxYyj1U+KyzzWD2vuKQ0vt1/lW/3X+VBSjoA7trSTPKtSSePfwUsnQ6OfQM7ZkJqgv7xOK3fAJ/J+pndhRAC/WjV1atXn/sZc0I8DwcHB7RabaaPNsrt97cEq3xWXILVI3GJaXxz4Crf/Stg1Smv+WcEy+XxL2fcTfhjClz457EFTnVg6BYoVS7fahVCFG46nU5OB4oCY2Fhke1IVW6/v+VUoHgu9rYWvP5CLYa3qqIfwTpwjXO34xmzOgSP8hpe861JJw8XVPYVYcBPcGY9bH0LNOXBtmxBly+EKETUavVzzbwuRGEgI1b5rLiNWP1XbGIq3+y/yncHrpHwzwhWXVf9CNYLHv+MYCXGQHqKPlwBJMXCrRD9tA1CCCFEIVRsH2kTExPDwIED0Wg0ODg4MGLECBISsn+USnJyMgEBAZQtWxY7Ozv8/f2feNp1eHg43bp1w9bWFmdnZ958803S09ON2qSkpPD2229TuXJlrKysqFKlCt9++63J97Eoc7C1ZEqn2ux7qz3j29eglKUZZyLiGf1DCN0X7ifwbBSKTZnHoQog8F1Y3RvWj4aH9wqueCGEEOI5FblgNXDgQM6cOUNgYCBbtmxh7969jB49Ots+kydPZvPmzaxbt449e/YQERFB7969DZ9nZGTQrVs3UlNTOXjwIN9//z0rV65kxowZRuvp27cvQUFBfPPNN4SFhfHjjz9Su3btPNnPoq5MKUve8KvN/qkdCGhf3RCwRq06Ro9F+9lxNko/X4iigGVpQAUnf4bFTeHkOv1yIYQQoogpUqcCz507h4eHB0ePHsXLywuAbdu20bVrV27evImrq+sTfeLi4nBycmLt2rX06dMHgPPnz1OnTh2Cg4Np0aIFf/75J927dyciIgIXFxcAli1bxtSpU7lz5w6WlpZs27aN/v37c+XKFRwdHXNdc0pKCikpKYb38fHxuLm5mfRUYFEQ8zCVr/ddYeXBaySmZgBQv4I9k3xr0sHdGdWtENg0AaLP6jvU7ATdvgQHtwKsWgghhNArlqcCg4ODcXBwMIQqAF9fX9RqNYcPH860T0hICGlpafj6+hqWubu7U6lSJYKDgw3rrV+/viFUAfj5+REfH8+ZM2cA2LRpE15eXnz22WdUqFCBWrVq8cYbb5CUlJRtzbNmzcLe3t7wcnMrmUHBsZQlb3V2Z//UDoxtVx1bSzNO3YpjxPfH6Ln4ADsT3FBG7Yb2b+unZLj4l35i0cs7C7p0IYQQIteKVLCKjIzE2dnZaJm5uTmOjo5ERkZm2cfS0hIHBwej5S4uLoY+kZGRRqHq0eePPgO4cuUK+/fv5/Tp0/z+++/MmzePX3/9lXHjxmVb8/Tp04mLizO8btwo2Y93cSxlydTO7ux7qz1j2uoD1smbcQxfeYyXlh9ll8tQlNF7wa05mFtD+UYFXbIQQgiRa4UiWE2bNg2VSpXt6/z58wVao06nQ6VSsWbNGpo1a0bXrl358ssv+f7777MdtbKyskKj0Ri9BJS1s2JaF33AerVtNWwszPj7ZhzDVh7lpV9j2NXyB5Th28H2n9OuigKnfoV0meNGCCFE4VUo5rGaMmUKQ4cOzbZNtWrV0Gq1REdHGy1PT08nJiYGrVabaT+tVktqaiqxsbFGo1ZRUVGGPlqtliNHjhj1e3TX4KM25cuXp0KFCtjb2xva1KlTB0VRuHnzJjVr1szVvgpjZe2smN6lDqNaV+OrvVdYFXydv2/EMuz7YzRyc2CSr4a2tZxQnfkdfhsBzh76x+JU9Mp55UIIIUQ+KxQjVk5OTri7u2f7srS0xNvbm9jYWEJCQgx9d+7ciU6no3nz5pmu29PTEwsLC4KCggzLwsLCCA8Px9vbGwBvb29OnTplFNoCAwPRaDR4eHgA0KpVKyIiIoymdrhw4QJqtZqKFSua9HiUROXsrJjetQ77prZndJtqWFuoCb0Ry9DvjtJ76UH2RFqg2JTVX9z+tS/8OQ1Ssp9mQwghhMhvRequQIAuXboQFRXFsmXLSEtLY9iwYXh5ebF27VoAbt26RceOHVm1ahXNmjUDYOzYsWzdupWVK1ei0WiYMGECAAcPHgT00y00atQIV1dXPvvsMyIjIxk0aBAjR47kk08+ASAhIYE6derQokUL3n//fe7evcvIkSNp27YtX331Va7rz4sJQoujOw9SWLH3Mj8cuk5ymv7ZYU0q2jGpVCCtry1EpQLsK0GPuVDDN/uVCSGEEM+pWN4VCLBmzRrc3d3p2LEjXbt2xcfHhxUrVhg+T0tLIywsjMTERMOyuXPn0r17d/z9/WnTpg1arZb169cbPjczM2PLli2YmZnh7e3NK6+8wuDBg/nggw8Mbezs7AgMDCQ2NhYvLy8GDhxIjx49WLBgQf7seAnjVNqKt7t5sPet9oz0qYqVuZrjNxMYHOZNH/sf2WfdASU2HFb7w473C7pcIYQQAiiCI1ZFnYxYPZvoB8ks33OF1Yeuk5KuH8Hy0sQxKXkJrYZ8iKp6+wKuUAghRHGW2+9vCVb5TILV84mOT2bZniusOfw4YDWtUoZJvrVoWb0sqot/gUtdsJfr3oQQQpiOBKtCSoKVaUTHJ7N0z2XWHA4n9Z+A1ayiLZNiPsTbLAzVCzPBawSoi9zZbiGEEIWQBKtCSoKVaUXFJ7N092XWHvlXwFKdY5L5b7SsUlo/NYOTPM9RCCHE85FgVUhJsMobkXHJLNtjHLCaq84yyWoT3u17QKtJYG5ZsEUKIYQosiRYFVISrPJWZFwyS3df4scj4aRm6H+1W6jPMMk5lBYTvgcziwKuUAghRFEkwaqQkmCVP27HJbF012V+OnKNVJ0KAO9qZZnkW5Pm1coWcHVCCCGKGglWhZQEq/wVEZvE0h1n+PlEtGEEq6WbNZMa6mjm80IBVyeEEKKokGBVSEmwKhgRsUks2X2Jn4/eIO2fgNVKc4dJvdvR1L1KwRYnhBCi0JNgVUhJsCpYt+7EsmTtr/xy24m0f55B7qPNYFLPVnhVlVOEQgghMifBqpCSYFU43DxziCXr/2Ldw4aGgNW6SmkmdamHZ2XHAq5OCCFEYSPBqpCSYFWIpKdyM3ARiw9EsC69NemPAlbNckzyrYVn5TIFXKAQQojCQoJVISXBqhCKPs+N395mSYwn6x42Il2n/79Em1pOTPKtSZNKErCEEKKkk2BVSEmwKqR0GZAYw43UUizedYlfQ24aAlbbfwJWYwlYQghRYkmwKqQkWBUN4Zs+YfHhGH7NaEsG+ucNtqvtxCTfWjRycyjY4oQQQuS73H5/yxNqhchEpaq1mG2/nl2Wr9PXbDdmKOwOu8NLiw8w7LsjhN6ILegShRBCFEIyYpXPZMSqCHl4D7ZPh5M/c13nzCKzV1if7MU/02DRwd2Z1zrWpKGMYAkhRLEnpwILKQlWRdDFHbBlEsTd4JrOhUXaj/j9ph0Z/1yD1dHdmdd8a9KgokOBlimEECLvSLAqpCRYFVEpCbDzQwhdC2MPci3dkYU7L/H7iZv8k6/wrePMax1rUb+ifcHWKoQQwuQkWBVSEqyKuMQYsH08gejVnd+yMMKdDWfu/ytguTDJtyb1KkjAEkKI4kIuXhciL/wrVHF5F1X3TubL8L7s8I2id2NX1CrYcS6K7gv3M2rVMU7fiiu4WoUQQuQ7CVZCPCuNK1RsBqkPqLZvMl8+nM6OoZXp1bgCahUEntUHrNGrjnEmQgKWEEKUBHIqMJ/JqcBiRpcBR7+GHe9D2kMws4S2b3G51igW7r7Kpr8jDKcI/eq68FrHWni4yn93IYQoauQaq0JKglUxFXsDtkyGS4H699Xaw+ANXIpOYOHOi2z6O4JH/0/rXFfLa741qVNe/vsLIURRIcGqkJJgVYwpCpz6FbZNha5fQL3eho8uRT9gQdAlNp98HLC61NMysaMELCGEKAokWBVSEqxKgOQ4sNKASqV/H/an/hRhjY5cin7A/KBLbPlXwOpaXx+w3LXy+yCEEIWVBKtCSoJVCfPwHixuCon3oOHL4Pcx2DpyMeoB84Mu8sep24aA1a1+eSZ2rEltbemCrVkIIcQTJFgVUhKsSphHE4seXg4oUMoJusyGur1BpeLCPwFr6z8BS6WCrvXL81rHmtRykYAlhBCFhQSrQkqCVQl14whsmgB3zuvf1+oC3eaAfQUAwiIfsOCfESzQB6xu/wSsmhKwhBCiwEmwKqQkWJVg6Smwfy7s/QJ0afrrsCaEgJ2zocn5yHgWBF1k66lIQB+wujdwZWKHGhKwhBCiAEmwKqQkWAmiz+lHr5zrwIsLM21yPjKe+Tsu8ufpxwGrRwNXJnasQQ1nCVhCCJHfiu0jbWJiYhg4cCAajQYHBwdGjBhBQkJCtn2Sk5MJCAigbNmy2NnZ4e/vT1RUlFGb8PBwunXrhq2tLc7Ozrz55pukp6cbtVmzZg0NGzbE1taW8uXLM3z4cO7du2fyfRTFnHMdGL4dOs9+vCw2HPbNgYw0ANy1Gpa+4smfr7Wmc10tigKb/o7ghbl7ee2nE1yKzv53XgghRMEocsFq4MCBnDlzhsDAQLZs2cLevXsZPXp0tn0mT57M5s2bWbduHXv27CEiIoLevR/PMZSRkUG3bt1ITU3l4MGDfP/996xcuZIZM2YY2hw4cIDBgwczYsQIzpw5w7p16zhy5AijRo3Ks30VxZjaDCxt9T8rCmyeBEEfwPK2cCvE0KxOeQ3LBnmydWJr/Oq6oCiwMTSCTnP3MOmnE1y+IwFLCCEKkyJ1KvDcuXN4eHhw9OhRvLy8ANi2bRtdu3bl5s2buLq6PtEnLi4OJycn1q5dS58+fQA4f/48derUITg4mBYtWvDnn3/SvXt3IiIicHFxAWDZsmVMnTqVO3fuYGlpyRdffMHSpUu5fPmyYd0LFy5k9uzZ3Lx5M9f7IKcCxRMUBU6tgz+nQlIMqNTQfCx0eBssSxk1PRMRx/wdF/nrrH7EVa2Cno0qMKFDDao52RVE9UIIUSIUy1OBwcHBODg4GEIVgK+vL2q1msOHD2faJyQkhLS0NHx9fQ3L3N3dqVSpEsHBwYb11q9f3xCqAPz8/IiPj+fMmTMAeHt7c+PGDbZu3YqiKERFRfHrr7/StWvXbGtOSUkhPj7e6CWEEZUKGvSF8Uehfl9QdHBoMSzxhss7jZrWdbVnxWAvtkzw4QUPF3QK/H7iFr5f7uH1n0O5IiNYQghRoIpUsIqMjMTZ2dlombm5OY6OjkRGRmbZx9LSEgcHB6PlLi4uhj6RkZFGoerR548+A2jVqhVr1qyhX79+WFpaotVqsbe3Z/HixdnWPGvWLOzt7Q0vNze3XO+vKGFKlQP/r+DldaCpCLHX4YdecHbjE03rVbDnq38Clm8dfcBa/yhg/RLK1bsPC2AHhBBCFIpgNW3aNFQqVbav8+fPF2iNZ8+e5bXXXmPGjBmEhISwbds2rl27xpgxY7LtN336dOLi4gyvGzdu5FPFosiq1QkCDkGzV6Fcbajpl2XTehXs+XqIF5vH++Bbx1kfsI7rA9aUX/7mmgQsIYTIV+YFXQDAlClTGDp0aLZtqlWrhlarJTo62mh5eno6MTExaLXaTPtptVpSU1OJjY01GrWKiooy9NFqtRw5csSo36O7Bh+1mTVrFq1ateLNN98EoEGDBpQqVYrWrVvz0UcfUb58+Uy3b2VlhZWVVbb7JsQTrEpD188gLQksrPXLMtIhaKb++qt/JhZ9pH5Fe74e0pRTN+OYt+MCQeej+e34TTaE3qJXY/01WJXLlnpyO0IIIUyqUIxYOTk54e7unu3L0tISb29vYmNjCQl5fNfUzp070el0NG/ePNN1e3p6YmFhQVBQkGFZWFgY4eHheHt7A/rrp06dOmUU2gIDA9FoNHh4eACQmJiIWm18uMzMzAAoQtf/i6LGwubxz4eXwsGFsLg5HP0GdLonmtevaM83Q5uyaXwrOrg7k6FT+DXkJh3m7OHNdX9z/Z6MYAkhRF4qUncFAnTp0oWoqCiWLVtGWloaw4YNw8vLi7Vr1wJw69YtOnbsyKpVq2jWrBkAY8eOZevWraxcuRKNRsOECRMAOHjwIKCfbqFRo0a4urry2WefERkZyaBBgxg5ciSffPIJACtXrmTUqFEsWLAAPz8/bt++zaRJk7K9cD4zclegeGaPJha9eVT/vlJLeHEBlKuZZZe/b8Qyb8cFdoXdAcBMrcK/SQXGt69JpbK2+VG1EEIUC8V25vWYmBjGjx/P5s2bUavV+Pv7s2DBAuzs9LeaX7t2japVq7Jr1y7atWsH6CcInTJlCj/++CMpKSn4+fmxZMkSo9OH169fZ+zYsezevZtSpUoxZMgQPv30U8zNH58tXbhwIcuWLePq1as4ODjQoUMHZs+eTYUKxqdlsiPBSjwXXQYc+Uo/51XaQzCzgrZvQavXwMwiy26h/wSs3f8ELHO1Cv8mFRnfoQZujhKwhBAiJ8U2WBV1EqyEScSGw5bJcGmH/n2jV+Cl7O9QBTgRfp95Oy6y58LjgNXHsyIB7SVgCSFEdiRYFVISrITJKAqc/AUCZ8DgDfpH5eTS8X8C1t5/Baz/eVVkXDsJWEIIkRkJVoWUBCthcukpYP6vO08PLgKXulC9fY5dQ67fZ96OC+y7eBd4FLDcCGhfnYplJGAJIcQjEqwKKQlWIk9FhMJX7fWztzd6BTp9CLaOOXYLuR7DvB0XDQHLwuxRwKpBBQebHHoLIUTxJ8GqkJJgJfJUygMI+hCOrAAUKOWsnw/L4yX9o3NycOyaPmDtv/Q4YPX1cmOcBCwhRAknwaqQkmAl8kX4Yf3UDHfD9O9rd4NuX4DmyQeVZ+botRjm7bjAgUv3AH3A6tfUjXHtauAqAUsIUQJJsCqkJFiJfJOeAvu+hH1zQJcGmgowMRTMLXO9iiNX9QHr4GV9wLI0U+sDVvvqlLeXgCWEKDkkWBVSEqxEvns0sWjD/tB05DOt4vCVe8zbcZHgK48DVv9mboxtJwFLCFEySLAqpCRYiQKhywBU8OixTJd3wq3jOU4s+l+Hrtxj3o4LHLoSA+gD1oBmboxtVwOtvXUeFC6EEIWDBKtCSoKVKHBpSbC4mX6SUZd68OJCqNDkqVYRfPkec3dc4MjVfwKWuZqXm1VibLvquGgkYAkhih8JVoWUBCtR4B5NLLptGiTFgEoNLcZB+7fB8unmrjp4+S7zAi9y5JoELCFE8SbBqpCSYCUKjYd39eHq1Dr9+zJVoMd8qNbuqVajKIphBOvotfsAWJmrebl5Jca2rY6zBCwhRDEgwaqQkmAlCp0L22HL6xB/U/9+1K6nPjUI+oB18PI95gZe4Nh1CVhCiOJFglUhJcFKFEopDyDoA3hwG/r+kKvJRLOiKAoHLulHsEL+FbAGNq/MmHbVcC4tAUsIUfRIsCqkJFiJQk2XAWoz/c8P78Jf70DHGbmeWPTfFEVh/6W7zA28wPHwWEAfsF5pUZlX20rAEkIULRKsCikJVqLIWP8qnPwJrDTwwvvQZOjj6RqegqIo7Lt4l7k7LnDin4BlbaHmleaVebVtdZxKW2W/AiGEKAQkWBVSEqxEkRF1Vj+x6K1j+veVW0GPBVCuxjOtTlEU9l7Uj2CF3ogF9AFrUAt9wCpnJwFLCFF4SbAqpCRYiSJFl6F/oHPQB5CWCGZW0G4qtJz4VBOL/puiKOy5cId5Oy4aApaNhRmDvCszuk21IhuwElPT8ZixHYCzH/hha2lewBUJIUwpt9/fTz+uL4QoOdRm0GIsjDsE1TtARoo+ZO378plXqVKpaFfbmd/HteS7YU1p6OZAUloGK/ZeofXsXczaeo57CSkm3AkhhMg/EqyEEDkrUxleWQ+9lutna28x5rlXqVKpaF/bmQ3jWvLd0KY0rGhPUloGy/dewWf2Lmb9KQFLCFH0SLASQuSOSqV/kPOr+8DaXr9MUfTXYV3Z/RyrVdHe3ZkNAa34dqgXDR4FrD1XaP3ZLj798zwxD1NNsw9CCJHHJFgJIZ7Ov+8MPPUrHF8Fq3rChgBIuv/Mq1WpVHRwd2FjQCu+GeJF/Qr2JKZmsGzPZXxm72T2NglYQojCT4KVEOLZ1fKDpqMAFYSuhkXN4MwG/UjWM1KpVHSs48Km8a34erAX9SpoSEzNYOnuy7SevZPPtp3nvgQsIUQhJcFKCPHsrDXQ7QsYvg3K1YKH0bBuCPw0EOIjnmvVKpUKXw8XNo/34avBXtR11fAwNYMlu/UjWJ9vP09sogQsIUThIsFKCPH8KrWAMfuhzVugNoewP+CXwc81cvWISqXiBQ8XtkzwYcUgTzzK6wPW4l2X8Zm9iy+2h0nAEkIUGhKshBCmYW4FHd6GV/dCxabQ6ePneubgf6lUKjrV1fLHRB+WD/KkTnkNCSnpLNp1CZ/Zu5jzVxhxiWkm254QQjwLmSA0n8kEoaJEUBTjUHVoqX6C0eeYWPS/dDqFv85GMW/HBc5HPgCgtJU5w1pVYYRPNextTbOd3JIJQoUo3mSCUCFEwfl3qIq7CTtm6icWXdEebh03ySbUahWd62nZOrE1y15pgru2NA9S0lmw8xI+s3fyZeAF4pJkBEsIkb8kWAkh8pamAvSYDzZlIOoUfN0R/noHUhNNsnp9wCrP1omtWTrwXwEr6CI+s3cyVwKWECIfSbASQuStRxOLBhyFen1A0cHBhbDU+7kmFv0vtVpFl/r6gLVkYBNqu5TmQXI68/8JWPN2XCA+WQKWECJvSbASQuQPOyfo8w0M+Fk/inX/Gqz5HzyINOlm1GoVXeuX58/XWrP45SbUcrHjQXI683ZcxOfTnczfcVEClhAiz8jVlUKI/FW7M1RuCUHvg50WSmvzZDNqtYpuDcrTpZ6WP09HMj/oAheiEpi74wLf7L/CyNbVGNaqCqWt8/cidyFE8VbkRqxiYmIYOHAgGo0GBwcHRowYQUJCQrZ9kpOTCQgIoGzZstjZ2eHv709UVJRRm4kTJ+Lp6YmVlRWNGjXKdD0nT56kdevWWFtb4+bmxmeffWaq3RKiZLHWQLc50OaNx8tuhcDPgyD+tkk39ShgbXutDYtebkxNZzvik9P5MvACPrN3sTDoIg9kBEsIYSJFLlgNHDiQM2fOEBgYyJYtW9i7dy+jR4/Ots/kyZPZvHkz69atY8+ePURERNC7d+8n2g0fPpx+/fpluo74+Hg6depE5cqVCQkJ4fPPP2fmzJmsWLHCJPslRIn06O5BRYE/psC5TbC4OYSsBJ3OpJtSq1V0b+DKtkltWDigMTWc7YhLSmPOPwFr0U4JWEKI51ek5rE6d+4cHh4eHD16FC8vLwC2bdtG165duXnzJq6urk/0iYuLw8nJibVr19KnTx8Azp8/T506dQgODqZFixZG7WfOnMmGDRsIDQ01Wr506VLefvttIiMjsbS0BGDatGls2LCB8+fPZ1lzSkoKKSkphvfx8fG4ubnJPFZC/FfUGdg0QT9yBVDZB15cAGWr58nmMnQKf5y6zfwdF7h85yEADrYWjGpdjSEtq2Bn9XRXSsg8VkIUb8VyHqvg4GAcHBwMoQrA19cXtVrN4cOHM+0TEhJCWloavr6+hmXu7u5UqlSJ4ODgp9p2mzZtDKEKwM/Pj7CwMO7fv59lv1mzZmFvb294ubm55XqbQpQoLnVhRCD4zQILW7i+H5Z4w74vIcP0I0lmahUvNnTlr8ltmd+/EdWdShGbmMbn28Pwmb2TxbsukZCSbvLtCiGKtyIVrCIjI3F2djZaZm5ujqOjI5GRmd9Z9GiEycHBwWi5i4tLln2yWo+Li8sT63j0WVamT59OXFyc4XXjxo1cb1OIEkdtBt7jYFwwVGsPGSn6i9zPbsyzTZqpVfRsVMEQsKr9K2C1nr2TJbsv8VAClhAilwpFsJo2bRoqlSrbV3an2wozKysrNBqN0UsIkYMyVWDQ7/DSMnDvDnWfvCbS1B4FrMDJbZnXrxHVypXifmIan23Tj2At3X1ZApYQIkeF4iKAKVOmMHTo0GzbVKtWDa1WS3R0tNHy9PR0YmJi0Gozv2Vbq9WSmppKbGys0ahVVFRUln2yWs9/7yR89P5p1iOEyCWVChoN0L8eSUmAn16G1lOgWts82ayZWsVLjSvQvUF5Np+MYEHQJa7efcjsbef5at8VRrepxmDvynINlRAiU4XiL4OTkxNOTk45tvP29iY2NpaQkBA8PT0B2LlzJzqdjubNm2fax9PTEwsLC4KCgvD39wcgLCyM8PBwvL29c12jt7c3b7/9NmlpaVhY6Oe9CQwMpHbt2pQpUybX6xFCPIcD8+DqHv2r8SDo9KH+UTl5wNxMTa/GFenRwJVNf0ewIOgi1+4l8umf5/lqrz5gDZKAJYT4j0JxKjC36tSpQ+fOnRk1ahRHjhzhwIEDjB8/nv79+xvuCLx16xbu7u4cOXIEAHt7e0aMGMHrr7/Orl27CAkJYdiwYXh7exvdEXjp0iVCQ0OJjIwkKSmJ0NBQQkNDSU1NBeDll1/G0tKSESNGcObMGX7++Wfmz5/P66+/nv8HQoiSquVEaDpK//OJH/RTM+Th9VegD1i9m1Rkx+ttmfO/hlQpa8u9h6nM+vM8rWfvYsXeyySmyilCIYRekZpuAfQThI4fP57NmzejVqvx9/dnwYIF2NnZAXDt2jWqVq3Krl27aNeuHaCfIHTKlCn8+OOPpKSk4Ofnx5IlS4xO4bVr1449e/Y8sb2rV69SpUoVQD9BaEBAAEePHqVcuXJMmDCBqVOnPlX9ub1dUwiRjfBD+qkZ7l7Qv3fvDl2/AE35PN90eoaODaERLNx5kev39A+SLmdnybBWVfl8exgg0y0IURzl9vu7yAWrok6ClRAmkpYM++bA/i9Blw4NB0CvZfm2+fQMHb+fuMXCnZcIj0k0+kyClRDFT7Gcx0oIIQwsrKHD2/DqXqjeAXzfz9fNm5up+Z+XG0FT2vJZnwZULGNj+OyzbWGkZ5h25nghRNEgI1b5TEashMgHGwL0M7a3nABm+fOQ5bikVBq+H2h436KaIwsHNMGptFW+bF8IkbdkxEoIUTKFH4LQ1fqJRb9qDxGh+bJZC7PHf05tLc04dCWGHgv3czw86yczCCGKHwlWQojixa25fmJRmzIQeQq+6gB/vQupiTn3NZFfXm1BdadSRMYn0295MD8cuo6cHBCiZJBgJYQoXh5NLBpwRD9ju5IBBxfA0pZwdW++lFDNyY6N433oUk9LWobCuxtO88a6kySnZeTL9oUQBUeClRCieLJzhv99BwN+gtKucP8q/D4W0lPyZ/NW5iwZ2IT/6+qOWgW/Hb9J7yUHuRGTfyNnQoj8J8FKCFG81e4CAYeh6UjoNgfM/7mYPB9OzalUKka3qc7qkc0pW8qSs7fj6b5wP7vDonPuLIQokiRYCSGKP2uNPlTV7vx42YnV8NNAiL+d55tvWb0cmyf40NDNgbikNIatPMqCoIvodHLdlRDFjQQrIUTJk5YMO2bC+S36x+KErARd3s475epgwy+vtuDl5pVQFPgy8AKjVh0jLiktT7crhMhfEqyEECWPhTUM3giuTSAlDja/BqtehHuX83SzVuZmfNKrPp/1aYCluZqg89G8uGg/527H5+l2hRD5R4KVEKJk0taDkTug08dgbgPX9unvHNw/FzLydhSpr5cb68e2pIKDDdfvJdJryQE2nLiVp9sUQuQPCVZCiJJLbQYtx8O4YKjWDtKTYcf7EHU6zzddr4I9Wyb40KaWE8lpOib9HMrMTWdITZdH4QhRlEmwEkIIx6owaAO8tBTavAmujR9/lod3D5YpZcl3Q5syoUMNAFYevMbLXx0iOj45z7YphMhbEqyEEAL+mVj0Zf2DnR+5eynPJxY1U6uY0qk2Xw32orSVOceu36fbwv0cuRqTZ9sUQuQdCVZCCJGVXR9D9Fn4vgdsmgBJsXm2qRc8XNg0wYfaLqW58yCFl786xLf7r8qjcIQoYiRYCSFEVnrMA68R+p+Pr4LFzeDspjzbXNVypfg9oCUvNnQlXafwwZazTPo5lMTU9DzbphDCtJ46WCUkJORFHUIIUfhY20P3L2HYn1C2BiREwS+D4OdX4EFknmzS1tKc+f0bMaO7B+ZqFRtDI+i1+CBX7z7Mk+0JIUzrqYOVvb09v/32W17UIoQQhVPlljDmALR+A9TmcG4z/P1jnm1OpVIx3Kcqa0e1oJydFWFRD3hx0X52nI3Ks20KIUzjqYOVoigsX76cVq1a4ePjw6RJkzh69Ghe1CaEEIWHhTV0fBdG74aGA8B7/OPP8mjW9mZVHfljog9elcvwIDmdkauOMeevMDLkUThCFFrPdI3ViRMnaNKkCT4+Ppw5c4bWrVvzxhtvmLo2IYQofLT1odcyMLPQv09Pha87Yh68ADMyTL45F401a0e1YGjLKgAs3HmJYSuPcv9hqsm3JYR4fubP0mnt2rW88MILhvcnT56kZ8+eVKhQgcmTJ5usOCGEKPRO/wYRx7GMOM5GyypMTRtt8k1YmquZ+WJdGrrZM339KfZeuEOPRftZ9oon9SrYm3x7Qohn99QjVo6Ojri5uRkta9CgAYsWLWLp0qUmK0wIIYqEhv2h5xIUawfqqa+x0fIdLHbOhLQkk2+qV+OK/D6uFZXL2nLzfhL+Sw+y7tgNk29HCPHsnjpYNWrUiO++++6J5TVq1CA8PNwkRQkhRJGhUkHjgSSNDmZLRgvMVTosDi38Z2LRfSbfXJ3yGjYF+NDR3ZmUdB1v/nqS//v9FCnppj8NKYR4ek8drD766CMWLFjAoEGDCA4O5uHDh0RHR/PJJ59QtWrVvKhRCCEKPztnxqdNZGTqFHR2Woi5Antm58kjcextLfhqsBeTfWuhUsHaw+H0W36I23GmHyUTQjydpw5WLVq04NChQ9y4cYPWrVuj0WgoX748v/76K3PmzMmLGoUQosjYofMkeXQwNB0FPebrR7QAMkw7yadareI135p8O7Qp9jYWhN6IpfuC/Ry8fNek2xFCPB2V8hzPS4iOjiYkJASdTkfz5s0pV66cKWsrluLj47G3tycuLg6NRlPQ5QghTCQxNR2PGdsBOPuBH7aW/7k36I83ICESun4BpbUm3Xb4vUTGrA7h7O141CqY1sWdUa2roXoU6oQQzy2339/PFazE05NgJUQJFHcLFjSCjFT9bO6dPoLGgx6PZplAUmoGb284xfrjtwDoWl/LZ30aYmf1TDd/CyH+I7ff3/KsQCGEyGv2FWDUTijfCJLj9A90/r4H3Ltssk3YWJox538N+fCleliYqdh6KpKXFh/gUrQ8hkyI/CTBSggh8oO2PowM0o9WmdvAtX36Owf3zzPZ9VcqlYpBLSrz02hvXDRWXIpO4KXFB9h2+rZJ1i+EyJkEKyGEyC9m5tByAow7CFXbQnoy7J8LSTEm3Yxn5TJsmdCa5lUdSUhJZ8zq48z68xzpGXnz6B0hxGPPHKxu3ryJ7l/Px9LpdPkyj1VMTAwDBw5Eo9Hg4ODAiBEjSEjIfqg7OTmZgIAAypYti52dHf7+/kRFGT/MdOLEiXh6emJlZUWjRo2eWMfu3bvp2bMn5cuXp1SpUjRq1Ig1a9aYcteEECWFYzUYvBF6LoZuc8DO+fFn6aZ5VI1TaSvWjGzOqNb6aXCW77nC4G+PcC8hxSTrF0Jk7pmDlYeHB9euXTO8v3PnTr7MYzVw4EDOnDlDYGAgW7ZsYe/evYwenf0jJCZPnszmzZtZt24de/bsISIigt69ez/Rbvjw4fTr1y/TdRw8eJAGDRrw22+/cfLkSYYNG8bgwYPZsmWLSfZLCFHCqFTQ+BWo3+fxsnNbYElzk00sam6m5u1uHix6uTG2lmYcvHyP7gv3E3oj1iTrF0I86ZnvCixdujR///031apVAyAqKory5csbjWKZ2rlz5/Dw8ODo0aN4eXkBsG3bNrp27crNmzdxdXV9ok9cXBxOTk6sXbuWPn30f8DOnz9PnTp1CA4OpkWLFkbtZ86cyYYNGwgNDc2xnm7duuHi4sK3336b632QuwKFEJlSFPi6I9wK0b9vMgRe+ABsHEyy+otRD3j1hxCu3H2IpZn+2YMDmrnJlAxC5FKB3BWY1/8HDQ4OxsHBwRCqAHx9fVGr1Rw+fDjTPiEhIaSlpeHr62tY5u7uTqVKlQgODn6ueuLi4nB0dMy2TUpKCvHx8UYvIYR4gkoFg34Hr+H698e/h8XN4dxmk6y+pktpNo5vhV9dF1IzdPzf76eY+ttJktPkUThCmNJTTXCyatUqw8/p6emsX78eZ2f9tQF+fn6mrSwTkZGRhu09Ym5ujqOjI5GRkVn2sbS0xMHBwWi5i4tLln1y45dffuHo0aMsX74823azZs3i/ffff+btCCFKEGt76D4X6vWBzRPh3iX4+RWo8+I/E4u6PNfqS1tbsOwVT5btucLn28/zy7GbnLv9gCUDm+DmaGuinRCiZHuqEavvvvvO8EpLS+PXX3/lu+++Y+XKlc9VxLRp01CpVNm+zp8//1zbMKVdu3YxbNgwvvrqK+rWrZtt2+nTpxMXF2d43bghT6IXQuSgSisYcwBaTwG1OZzbBLf/NsmqVSoVY9tVZ9Xw5pSxteDUrTh6LNrP3gt3TLJ+IUq6pxqx2rVrl+Hn0qVLs3btWqNrrJ7VlClTGDp0aLZtqlWrhlarJTo62mh5eno6MTExaLWZPyJCq9WSmppKbGys0ahVVFRUln2ys2fPHnr06MHcuXMZPHhwju2trKywsrJ66u0IIUo4C2voOAPq9oLzW6FWp8efpSXrP38OPjXLsWVia8auDuHkzTiGfHeENzrVZmzb6qjVct2VEM+qUDzrwMnJCScnpxzbeXt7ExsbS0hICJ6engDs3LnT8KzCzHh6emJhYUFQUBD+/v4AhIWFER4ejre391PVuXv3brp3787s2bNzvBNRCCFMQltf/3ok/jasaAfe46BFgH5urGdUwcGGX171ZuamM/x09Aafbw8j9EYsc/o2RGNt8fy1C1ECmfTi9bx+7GCdOnXo3Lkzo0aN4siRIxw4cIDx48fTv39/wx2Bt27dwt3dnSNHjgBgb2/PiBEjeP3119m1axchISEMGzYMb29vozsCL126RGhoKJGRkSQlJREaGkpoaCipqfo5ZXbt2kW3bt2YOHEi/v7+REZGEhkZSUyMaSf2E0KIbB1fpX+Yc+AM+LoD3D75XKuztjDjU/8GfNq7PpbmagLPRtFz0QHCIh+YqGAhShjlGX3yySfK/fv3De8TEhKUmTNnPuvqcu3evXvKgAEDFDs7O0Wj0SjDhg1THjx4YPj86tWrCqDs2rXLsCwpKUkZN26cUqZMGcXW1lbp1auXcvv2baP1tm3bVgGeeF29elVRFEUZMmRIpp+3bdv2qeqPi4tTACUuLu5ZD4EQoiTT6RTl+A+KMstNUd7TKMrMMooS+J6ipCY+96r/vnFfaTkrSKk8dYvi/s6fyqbQW89frxDFRG6/v595HivxbGQeKyGESTyIgj/fgrMb9O8dq8OLC6CKz3OtNuZhKhN/PMH+S3cBGOFTlWld3LEwkyegiZKtQOaxEkIIkU9Ku0Df76H/WihdHmIuw/k/nnu1jqUs+X54M8a1qw7AN/uvMvDrw0Q/SH7udQtREsiIVT6TESshhMklx8H+edDmDbAspV+W+vDxz89o2+lI3lj3Nwkp6bhorFgysAmelbOfFFmI4kpGrIQQoqSwtgff9x4HKZ0OfugNvwzWnzJ8Rp3radk4vhU1ne2Iik+h3/JDfH/wWp7fqCREUSbBSgghipuI43DzKJzdCIub6u8kfMYwVN3Jjg0BrejWoDzpOoX3Np3h9V/+JilVHoUjRGae6VRgeHg4e/fuxcrKisaNG1OjRo28qK1YklOBQoh8cfskbJoAt0P176u2gR7zwbHaM61OURS+2X+VWX+eJ0On4K4tzfJBnlQu+3ynG4UoKnL7/f3UwWrBggW8/vrr2NraolKpSEhIwNPTk6+//poGDRo8d+HFnQQrIUS+yUiHQ0tg1yeQngTmNtB++nNNLHroyj3Grz3O3YRUNNbmzOvfiA7uz/cMQyGKgjy7xurDDz9k2rRpxMbGEhcXR1hYGD4+Pnh7e7N///7nKloIIYQJmZlDq4kw7qB+xCo9CU6tQz8N37NpUa0sWya0pnElB+KT0xm+8hhzAy+g08l1V0LAM4xY2djYcOrUqSdO/3300Uds3ryZw4cPm7TA4kZGrIQQBUJRIHQNuNQF18b6ZempoGSAhc1Try41XceHW87yw6HrALSv7cTcfo1wsLU0ZdVCFBp5NmLVoEEDgoODn1jet29fTp58vkcrCCGEyCMqFTR+5XGoAtg3B5a2gmtPf7bB0lzNhy/VY87/GmJlrmZX2B16LNrPmYg4ExYtRNHz1MFqzpw5TJkyhZ9//tnoltvDhw9Ts2ZNkxYnhBAij6Qlw99r9ROLruwGm1/Tz4f1lPw9K7J+XEvcHG24EZNE7yUHWX/8Zh4ULETR8Ex3BW7dupUxY8aQnJxMo0aNSE1N5fTp0/zwww906dIlL+osNuRUoBCi0EiO0z/MOWSl/r2dFrrNgTrdn3pVsYmpTPo5lN1hdwAY1KIy73b3wNJcZvURxUOe3RX4SEpKCkFBQezevZsTJ05w4sQJ7t+/j6OjI/Xr16dBgwbMmzfvWesvtiRYCSEKnWv7YdNE/egVgEdP6PK5/rE5T0GnU5gXdJEFQRcBaFLJgSUDPdHaW5u6YiHyXZ4Hq8zcuHGD0NBQQ9D6/fffTbXqYkOClRCiUEpLgj2z4cACMLPU30n4jHNeBZ2LYvLPocQnp1POzpJFLzehRbWyJi5YiPxVIMFK5EyClRCiULt9Eu6EQYP/PV6WGAO2T/eMwOv3HvLqDyGcj3yAmVrF9C7ujPCpikqlMnHBQuQPeVagEEKIp1e+gXGoun4Q5taDgwv1E47mUuWypfh9XCteauRKhk7hoz/OMeHHEzxMyf06hCiKJFgJIYTI2smfIe0h/PUOfOMLkady3dXG0oy5/Rrx/ot1MVer2HLyNr2WHODKnYQ8LFiIgiXBSgghRNa6z4MXF4G1PUScgOVtYcf7+ukackGlUjGkZRV+Gt0C59JWXIhK4MVFB9h+JjJv6xaigEiwEkIIkTWVCpoMgoAjUOdF/Uzt+7+EZa3g+pOTRWfFq4ojWyb60KyKIwkp6bz6QwifbdM/0FmI4kSClRBCiJyV1kK/H6Dfav18V/cuQez1p1qFc2lr1oxqzvBWVQFYsvsyQ787QszD1LyoWIgCIXcF5jO5K1AIUeQlxULoWmgxVj+iBZBwB+yccr2KjaG3mPbbKZLSMqjgYMPSV5rQoKJDnpQrhCnIXYFCCCHyho0DeI97HKqSYmGZD/wyBBKic7WKno0q8HtAS6qUteVWbBJ9lgXz89HwPCtZiPwiwUoIIcTzubYPHt6BsxtgUVM4sRpycTLEXath0wQffOu4kJquY+pvp5i+/iTJaRl5X7MQeUSClRBCiOdTpweM3gXlG0JyLGwMgFU9IeZqjl011hasGOTJG51qoVLBj0du0Hd5MLdik/K+biHygAQrIYQQz698Qxi5E174AMyt4eoeWOKtn1g0h9ErtVrF+A41WTmsGQ62Fpy8GUePhfs5cOluPhUvhOlIsBJCCGEaZubQ6jUYexCqtoH0JIg8/fharBy0reXE5vE+1KugIeZhKoO+OczS3ZeRe6xEUSJ3BeYzuStQCFEiKIr+zsHaXR4/ZzAxBixs9K9sJKdl8O6G06wLuQmAX10XvvhfQ0pbW+R11UJkSe4KFEIIUXBUKmg88HGoUhT9tVfLfODagWy7WluY8VmfBnzSqz6WZmq2n4mi5+IDXIx6kA+FC/F8JFgJIYTIewlRcOu4fmLRlV1h8yRIjsuyuUql4uXmlfhljDfl7a25cuchPRcf4I+Tt/OvZiGegQQrIYQQea+0FgIOQ5Mh+vch38Hi5nB+a7bdGrk5sHmCD97VypKYmkHA2uN8svUc6Rm6fChaiKcnwUoIIUT+sHGAFxfAkC3gWA0e3IafBsC6ofrrr7JQzs6KH0Y049W21QBYsfcKr3xzmLsJKflTtxBPQYKVEEKI/FW1tf7OwVaTQGUGESf0UzRkw9xMzfQudVg6sAmlLM04dCWG7gv2czz8fv7ULEQuFblgFRMTw8CBA9FoNDg4ODBixAgSEhKy7ZOcnExAQABly5bFzs4Of39/oqKijNpMnDgRT09PrKysaNSoUbbru3TpEqVLl8bBweE590YIIUooCxt44X39xKK9VoClrX65LgNib2TZrUv98mwc34rqTqWIjE+m3/Jgfjh0XaZkEIVGkQtWAwcO5MyZMwQGBrJlyxb27t3L6NGjs+0zefJkNm/ezLp169izZw8RERH07t37iXbDhw+nX79+2a4rLS2NAQMG0Lp16+faDyGEEOgnFq3U/PH7o9/oH4tzcCFkpGfapYZzaTaO96FLPS1pGQrvbjjNG+vkUTiicChS81idO3cODw8Pjh49ipeXFwDbtm2ja9eu3Lx5E1dX1yf6xMXF4eTkxNq1a+nTpw8A58+fp06dOgQHB9OiRQuj9jNnzmTDhg2EhoZmWsPUqVOJiIigY8eOTJo0idjY2GxrTklJISXl8XUA8fHxuLm5yTxWQgjxX4oCP70MYf9c0O7aGF5cCNr6WTRXWLH3CrO3nUenQF1XDcte8cTN0TYfixYlRbGcxyo4OBgHBwdDqALw9fVFrVZz+PDhTPuEhISQlpaGr6+vYZm7uzuVKlUiODj4qba/c+dO1q1bx+LFi3PdZ9asWdjb2xtebm5uT7VNIYQoMVQq6L8WeiwAK3v9tVcr2kHQB5CWnElzFa+2rc7qEc0pW8qSMxHxdF+4n91h0flfuxD/KFLBKjIyEmdnZ6Nl5ubmODo6EhkZmWUfS0vLJ66HcnFxybJPZu7du8fQoUNZuXLlU400TZ8+nbi4OMPrxo2srx0QQogST6UCzyEw/oj+4c66dNg3B5a1gpshmXZpWaMcmyf40NDNgbikNIatPMqCoIvodEXmhIwoRgpFsJo2bRoqlSrb1/nz5wu0xlGjRvHyyy/Tpk2bp+pnZWWFRqMxegkhhMhBaS30Ww19fwA7F4i5CmZZP9LG1cGGX15twcvNK6Eo8GXgBUatOkZcUlo+Fi0EmBd0AQBTpkxh6NCh2bapVq0aWq2W6GjjId709HRiYmLQarWZ9tNqtaSmphIbG2s0ahUVFZVln8zs3LmTTZs28cUXXwD6c/s6nQ5zc3NWrFjB8OHDc70uIYQQueTxov6Bztf2QfkGj5ffvQTlahg1tTI345Ne9Wnk5sA7G04TdD6aFxftZ9krntQpL/+oFfmjUAQrJycnnJyccmzn7e1NbGwsISEheHp6AvrAo9PpaN68eaZ9PD09sbCwICgoCH9/fwDCwsIIDw/H29s71zUGBweTkfH4jpONGzcye/ZsDh48SIUKFXK9HiGEEE/JxkF/WvCRqDOwvI1+WZfPwM74EpG+Xm7U0WoYszqE6/cS6bXkALP9G9CzkfytFnmvUJwKzK06derQuXNnRo0axZEjRzhw4ADjx4+nf//+hjsCb926hbu7O0eOHAHA3t6eESNG8Prrr7Nr1y5CQkIYNmwY3t7eRncEXrp0idDQUCIjI0lKSiI0NJTQ0FBSU1MN265Xr57hVaFCBdRqNfXq1aNMmTL5fzCEEKKkunFYfwfhmd/1UzOcWKN//y/1K9qzZYIPrWuWIzlNx2s/hTJz0xnS5FE4Io8VqWAFsGbNGtzd3enYsSNdu3bFx8eHFStWGD5PS0sjLCyMxMREw7K5c+fSvXt3/P39adOmDVqtlvXr1xutd+TIkTRu3Jjly5dz4cIFGjduTOPGjYmIiMi3fRNCCJELXsNh1E79NAzJsbBxHPzwkv46rH8pU8qSlcOaMb69/pThyoPXGLDiENHxT95hKISpFKl5rIqD3M6DIYQQIgcZaRC8GHbPgvRkMLeBTh9Cs1FPNA08G8XrP4fyICUdp9JWLBnYhKZVHAugaFFUFct5rIQQQggDMwvwmaR/7mCV1pCeBBmpmTZ9wcOFTRN8qO1SmjsPUhiw4hDfHbgqj8IRJicjVvlMRqyEECIPKAqc3ai/oF1tpl927zJoKoDF4wc8J6amM/W3U2z+W3+ZR89GrszqXR9by0JxL5coxGTESgghRMmhUkHdlx6HqrQkWPM/WOYD1w8amtlamrOgfyNmdPfATK1iY2gEvZcc5NrdhwVTtyh2JFgJIYQofmKuQGoC3LsI33WBLa9DcjygfxTOcJ+q/DiqBeXsrDgf+YAei/az42xUARctigMJVkIIIYofl7oQcASaDNa/P/YNLG4OYX8amjSr6sgfE33wrFyGB8npjFx1jDl/hZEhj8IRz0GClRBCiOLJxgFeXAhDNkOZqvAgAn7sD+uG6k8VAi4aa34c1YKhLasAsHDnJYatPEpsYuYXwQuREwlWQgghireqbWBcMLR6DVRmkJIA5o8vaLc0VzPzxbrM7dcQaws1ey/cofvC/Zy+FVeARYuiSu4KzGdyV6AQQhSgiFCwLQsObvr3SfchOQ7KVAHgbEQ8Y1aHEB6TiJW5mo9eqsf/vNwKrFxReMhdgUIIIcR/uTZ6HKoAtr8DS7z1E43qMvBw1bB5vA8d3J1JSdfx5q8nefv3U6SkZ2S5SiH+TYKVEEKIkik9BWKvQ1oibP8/+NoXIk9jb2vB14O9mOxbC5UK1hwOp9/yQ9yOSyroikURIMFKCCFEyWRuBYM3QY/5YGUPEcdhRVsI+hB1Rgqv+dbk26FN0VibE3ojlu4L9nPw8t2CrloUchKshBBClFxqNXgOhYDD4N4ddOmw7wv9xKLR52hf25ktE1rjUV7DvYepvPL1YVbsvSyPwhFZkmAlhBBCaMpD/zXQdxXYuUByrP5/gUplbfltbEt6N6mAToFPtp4nYO1xElLSC7ZmUSjJXYH5TO4KFEKIQi7pvv45gxW99O8VBW4eRanYlNWHw/lg8xnSMhRqONux7BVPajjbFWy9Il/IXYFCCCHEs7Ap8zhUAZz+Db55AdWvwxlUz5afRnvjorHiUnQCLy0+wLbTtwuuVlHoSLASQgghshMbrp9Y9Mx6WNwUz/t/smW8D82rOpKQks6Y1ceZ9ec50jN0BV2pKAQkWAkhhBDZaf06jNoJ2vr604QbxuK0cQBr+mgZ6VMVgOV7rjD42yPcS0gp4GJFQZNgJYQQQuTEtRGM2gW+M/WPw7m8E/NlLXmn4t8serkxtpZmHLx8j+4L9xN6I7aAixUFSYKVEEIIkRtmFuAzGcYehCqt9ROL2legewNXNgS0olq5UtyOS6bvsmDWHg6XKRlKKLkrMJ/JXYFCCFEM6HQQHgxVWhkWxV86xBsH1Px17g4Afb0q8kHPelhbmBVUlcKE5K5AIYQQIq+o1UahivvX0fzci+XxAbzV3Aa1Cn45dpP/LQvm5v3EgqtT5DsJVkIIIcTzirsJlnao7l1g3N+9WFX3BGVszDl1K44eC/ez7+Kdgq5Q5BMJVkIIIcTzqtIKxh+BxoMA8Ln0OVts36dBObifmMbgb4+weNcldDq5+qa4k2AlhBBCmIJNGei5SP9g5zJVqPDwDL88GEJ/p+soCny+PYxXV4cQn5xW0JWKPCTBSgghhDClam1hbDC0nIi1OoNPG93l0971sTRTE3g2ip6LDhAW+aCgqxR5RO4KzGdyV6AQQpQgt/+GcrXBwpq/b8QybtVhbj1Ix8bCjM/6NKBHQ9eCrlDkktwVKIQQQhS08g3BwhqAhq52bHach4/ZWZLSMpjw4wk+3HKWNHkUTrEiwUoIIYTID8lxOFqr+N78Y8aabQTgm/1XGfj1YaIfJBdwccJUJFgJIYQQ+aFUWRiyBbMec5lq9yfLLL7EjiSOXI2hx4L9hFyPKegKhQlIsBJCCCHyi1oNXsMg4DCdPVzYaPkONVU3iXqQQv/lh1gVfE0ehVPEFblgFRMTw8CBA9FoNDg4ODBixAgSEhKy7ZOcnExAQABly5bFzs4Of39/oqKijNpMnDgRT09PrKysaNSoUabrURSFL774glq1amFlZUWFChX4+OOPTbVrQgghSgqNK/RfQ/X+s9lQZgHd7C6QplOYsfEMr//yN0mpGQVdoXhGRS5YDRw4kDNnzhAYGMiWLVvYu3cvo0ePzrbP5MmT2bx5M+vWrWPPnj1ERETQu3fvJ9oNHz6cfv36Zbme1157ja+//povvviC8+fPs2nTJpo1a/bc+ySEEKIEUqnAoyelJuxn0biXeLtrHczUKn4/cYve87Zz/d7Dgq5QPIMiNd3CuXPn8PDw4OjRo3h5eQGwbds2unbtys2bN3F1ffK21bi4OJycnFi7di19+vQB4Pz589SpU4fg4GBatGhh1H7mzJls2LCB0NDQJ7bdoEEDTp8+Te3atZ95H2S6BSGEEFkJvnyPCSv3cTfNAo1ZKvP/V4/2jWoWdFmCYjrdQnBwMA4ODoZQBeDr64tarebw4cOZ9gkJCSEtLQ1fX1/DMnd3dypVqkRwcHCut71582aqVavGli1bqFq1KlWqVGHkyJHExGR/sWFKSgrx8fFGLyGEECIz3tXLsrn5GRqrLhKfYcnwn84z94f16GRKhiKjSAWryMhInJ2djZaZm5vj6OhIZGRkln0sLS1xcHAwWu7i4pJln8xcuXKF69evs27dOlatWsXKlSsJCQkxjIJlZdasWdjb2xtebm5uud6mEEKIkqd897f5aWwbBtmFoKBm/hkrRsxaTtztKwVdmsiFQhGspk2bhkqlyvZ1/vz5Aq1Rp9ORkpLCqlWraN26Ne3ateObb75h165dhIWFZdlv+vTpxMXFGV43btzIx6qFEEIURVaVmvDh9OnMaRiBFansSqhEjwV7OHPgj4IuTeTAvKALAJgyZQpDhw7Ntk21atXQarVER0cbLU9PTycmJgatVptpP61WS2pqKrGxsUajVlFRUVn2yUz58uUxNzenVq1ahmV16tQBIDw8PMvrrqysrLCyssr1doQQQggAzCzwHzCK2g1OM/ank4SnOdP7TxWzbG7Su0nFgq5OZKFQBCsnJyecnJxybOft7U1sbCwhISF4enoCsHPnTnQ6Hc2bN8+0j6enJxYWFgQFBeHv7w9AWFgY4eHheHt757rGVq1akZ6ezuXLl6levToAFy5cAKBy5cq5Xo8QQgjxNOrVrcfmqTV4bfVB9lxL4vVf/ib0RizvVLuCZR0/MJd/vBcmRequQIAuXboQFRXFsmXLSEtLY9iwYXh5ebF27VoAbt26RceOHVm1apVhKoSxY8eydetWVq5ciUajYcKECQAcPHjQsN5Lly6RkJDAsmXL2LVrFz///DMAHh4eWFpaotPpaNq0KXZ2dsybNw+dTkdAQAAajYa//vor1/XLXYFCCCGeRYZOYX7QRRYEXQSgieoCS7Sb0fb6BCq1yKG3eF65/f4ucsEqJiaG8ePHs3nzZtRqNf7+/ixYsAA7OzsArl27RtWqVdm1axft2rUD9BOETpkyhR9//JGUlBT8/PxYsmSJ0anAdu3asWfPnie2d/XqVapUqQJAREQEEyZM4K+//qJUqVJ06dKFOXPm4OjomOv6JVgJIYR4HkHnopj04zEepEI5YllkuZAWzX3A9z2wKl3Q5RVbxTZYFXUSrIQQQjyva3cfMuaHI5yPSsSMDKabr2VEmZOoesyFWn4FXV6xVCznsRJCCCEEVClXit8D2vBSI1cyMOOj9EFMuNebh2sGwZ/TCrq8Ek2ClRBCCFEE2ViaMbdfI2b28MBcrWKLzpteqR9wpUzub8wSpifBSgghhCiiVCoVQ1tV5afRLXAubcUFxY2ef1rz15l/JsC+sgfuXy/YIksYCVZCCCFEEedVxZEtE3xoVsWRBynpjP4hhM82nyDjl2GwxBsOLQVdRkGXWSJIsBJCCCGKAWeNNWtGNWdYqyoALDkQwdDUN4lJVcO2afBNJ4g+V7BFlgASrIQQQohiwsJMzXs96jK/fyNsLMzYl1iJHmZLOGVeD24dg2WtYdcnkJ5S0KUWWxKshBBCiGKmZ6MK/B7QkiplbbmVaIZ/0tv8XG4i6NJgz2xY3gZSEwu6zGJJgpUQQghRDLlrNWwc74NvHWdSMxSm3mzB9Io/kGKrhUreYGlb0CUWSxKshBBCiGLK3saCFYO8eKNTLVQq+PGSGX2tl3Or6f89bnT/GlzYXmA1FjcSrIQQQohiTK1WMb5DTVYOa4aDrQV/Rzykx1d/c+DSXVAU2PwarO0Lv46Ah3cLutwiT4KVEEIIUQK0reXE5vE+1HXVEPMwlUHfHGbp7osozvVApYbTv8KipvD3T/rAJZ6JBCshhBCihHBztOW3sS35n2dFdArM3n6RMXd68WBwILjUg6QY+P1VWNMHYsMLutwiSYKVEEIIUYJYW5jxWZ8GfNyrHhZmKrafiaLn+gQuvrQFOs4AMyu4tAMWt4AbRwu63CJHgpUQQghRwqhUKgY2r8wvr3pT3t6aK3ce0nPpYf6wfxnGHoBKLcG+ApRvUNClFjkSrIQQQogSqnGlMmye4IN3tbIkpmYQsPY4nxxJJ33QZhiyGcyt9A0z0uHIVzKxaC5IsBJCCCFKsHJ2Vvwwohmvtq0GwIq9V3jlu6PcVZV53Ch4EWx9Qz9ze/jhAqq0aJBgJYQQQpRw5mZqpnepw9KBTShlacahKzF0X7Cf4+H39Q0cq0EpJ7gbBt/6wdY3IeVBwRZdSEmwEkIIIQQAXeqXZ+P4VlRzKkVkfDL9lgfzw6HrKHV6QMARaDQQUODICv3F7Rf+KuiSCx0JVkIIIYQwqOFcmo0BrehST0tahsK7G07zxrqTJFvYw0tLYNAGcKgM8Tdh7f9g7xcFXXKhIsFKCCGEEEZKW1uwZGATpnVxR62C347fxH/pQW7EJEL19jAuGLzHg7kNuHcv6HILFZWiyPSq+Sk+Ph57e3vi4uLQaDQFXY4QQgiRrYOX7jL+xxPEPEzF3saCBQMa07aWk/7DhDtg5/S48YnVULUNOFQqmGLzUG6/v2XESgghhBBZalmjHFsm+NDQzYG4pDSGfneEhUEX0ekU41B1KwQ2TdBfe3VoGegyCq7oAiTBSgghhBDZcnWw4ZdXW/By80ooCswJvMDoH44Rl5T2uJGVPbg1h7SHsG2q/u7B6HMFV3QBkWAlhBBCiBxZmZvxSa/6fNanAZbmanaci6bnov2cj4zXNyhXA4ZuhW5zwLI03Dyqn/dq1yclamJRCVZCCCGEyLW+Xm78NqYlFRxsuHYvkZcWH2Bj6C39h2o1NB0JAYehVhfQpcGe2bCyG5SQS7olWAkhhBDiqdSvaM+WCT60rlmO5DQdr/0UysxNZ0jL0Okb2FeAAT9Cn+/0E4s2ehlUqoItOp/IXYH5TO4KFEIIUVxk6BTmBl5g0a5LADStUobFLzfBWWP9uFHSff31V+p/xnKu7IGMVKj5QgFU/OzkrkAhhBBC5CkztYo3/GqzYpAnpa3MOXrtPt0W7ufotZjHjWzKPA5VKQmwMQDW9IHfRsHDuwVTeB6SYCWEEEKI59KprpaN41tRy8WOOw9SGLDiEN8duMoTJ8VUavDoqf/fU7/A4mZw8heTXH+VmJpOlWl/UGXaHySmpj/3+p6VBCshhBBCPLdqTnZsCGhFj4aupOsU3t98lkk/hxqHHEtb8PsYRu4Al3qQeA/Wj4I1/4PY8IIr3oSKXLCKiYlh4MCBaDQaHBwcGDFiBAkJCdn2SU5OJiAggLJly2JnZ4e/vz9RUVFGbSZOnIinpydWVlY0atQo0/Vs376dFi1aULp0aZycnPD39+fatWsm2jMhhBCiaLO1NGdB/0a8290DM7WKjaER9F5ykGt3Hxo3rOAJo3dDh3fAzBIuBeonFr1/vUDqNqUiF6wGDhzImTNnCAwMZMuWLezdu5fRo0dn22fy5Mls3ryZdevWsWfPHiIiIujdu/cT7YYPH06/fv0yXcfVq1fp2bMnHTp0IDQ0lO3bt3P37t1M1yOEEEKUVCqVihE+VVk7sjnl7Kw4H/mAHov2s+Os8YAGZhbQ5k0YcwAqeUONjlCmcsEUbUJF6q7Ac+fO4eHhwdGjR/Hy8gJg27ZtdO3alZs3b+Lq6vpEn7i4OJycnFi7di19+vQB4Pz589SpU4fg4GBatGhh1H7mzJls2LCB0NBQo+W//vorAwYMICUlBfU/F+Ft3ryZnj17kpKSgoWFRa72Qe4KFEIIUVJExSczbs1xQq7fB2BChxpM8q2Fmfo/Uy/odPoZ261K698n3IHj30PLCWBulattJaam4zFjOwBnP/DD1tLcZPsBxfSuwODgYBwcHAyhCsDX1xe1Ws3hw4cz7RMSEkJaWhq+vr6GZe7u7lSqVIng4OBcb9vT0xO1Ws13331HRkYGcXFx/PDDD/j6+mYbqlJSUoiPjzd6CSGEECWBi8aaH0e1YIi3fiRq4c5LDF95lNjEVOOGavXjUAX6R+Ls/BCWt4EbR/Kx4udXpIJVZGQkzs7ORsvMzc1xdHQkMjIyyz6WlpY4ODgYLXdxccmyT2aqVq3KX3/9xf/93/9hZWWFg4MDN2/e5Jdffsm236xZs7C3tze83Nzccr1NIYQQoqizNFfzfs96zO3XEGsLNXsu3KH7wv2cvhWXdSf3bmBbDu6ch286wda39FM1FAGFIlhNmzYNlUqV7ev8+fMFWmNkZCSjRo1iyJAhHD16lD179mBpaUmfPn2evJ30X6ZPn05cXJzhdePGjXysWgghhCgcejWuyPqxrajkaMvN+0n4Lz3IryE3M29czx/GH4WGLwMKHFkOS1rAxcB8rflZmPYE5DOaMmUKQ4cOzbZNtWrV0Gq1REdHGy1PT08nJiYGrVabaT+tVktqaiqxsbFGo1ZRUVFZ9snM4sWLsbe357PPPjMsW716NW5ubhw+fPiJa7UesbKywsoqd+eHhRBCiOLMw1XD5vE+TP4llJ3no3lj3d+cCL/PjB4eWJmbGTe2dYReS6F+H9gyST8dw5o+0HMJNB5YIPXnRqEIVk5OTjg5OeXYztvbm9jYWEJCQvD09ARg586d6HQ6mjdvnmkfT09PLCwsCAoKwt/fH4CwsDDCw8Px9vbOdY2JiYmGi9YfMTPT/xLodLpcr0cIIYQoyextLfh6sBcLd15iXtAF1hwO50xEPEtfaUJ5e5snO9ToCOMOwc6P4dwmqNMj/4t+CoXiVGBu1alTh86dOzNq1CiOHDnCgQMHGD9+PP379zfcEXjr1i3c3d05ckR/sZu9vT0jRozg9ddfZ9euXYSEhDBs2DC8vb2NRpkuXbpEaGgokZGRJCUlERoaSmhoKKmp+gvsunXrxtGjR/nggw+4ePEix48fZ9iwYVSuXJnGjRvn/8EQQgghiii1WsVrvjX5dkhTNNbmhN6IpfuC/Ry8nMUjbixLQedPIOAwWP9zR56iwM6PILZwXWJTpIIVwJo1a3B3d6djx4507doVHx8fVqxYYfg8LS2NsLAwEhMTDcvmzp1L9+7d8ff3p02bNmi1WtavX2+03pEjR9K4cWOWL1/OhQsXaNy4MY0bNyYiIgKADh06sHbtWjZs2EDjxo3p3LkzVlZWbNu2DRubTBK2EEIIIbLV3t2ZLRNaU6e8hnsPUxn0zRFW7L2c9bXLlqUe/3xiNez9XH/t1eEVoMvIn6JzUKTmsSoOZB4rIYQQwlhSagZv/36K9SduAdCtfnlm92mAnVU2VyzduQCbJ0K4fuqkjApN8bvyPy4pFWUeKyGEEEKUXDaWZszp25APe9bFwkzFH6du89LiA1yKzmaKBadaMHQrdJsDlqUxu3WUrZbTaav+O/8Kz4QEKyGEEEIUOJVKxSDvKvw02hsXjRWXohN4afEBtp2+nXUntRqajoSAQ2TU8COaMhzV1c6/ojMrqUC3LoQQQgjxL56Vy7BlQmuaV3UkISWdMauP8+mf50nPyOYOfPuKpPxvDb1SPiAR6/wrNhMSrIQQQghRqDiVtmL1yOaM9KkKwLI9lxn87RHuJaRk3Uml4g4O+VNgNiRYCSGEEKLQsTBT8053DxYOaIytpRkHL9+jx8L9hN6ILejSsiXBSgghhBCFVo+GrmwIaEW1cqWIiEum77JgfjwSXtBlZUmClRBCCCEKtVoupdkwvhWdPFxIzdAxff0ppv56kuS0wjF31b9JsBJCCCFEoaextmDZK5686VcbtQp+PnaD/y0L5ub9xJw75yMJVkIIIYQoEtRqFQHta/D98GaUsbXg1K04eizcz76Ldwq6NAMJVkIIIYQoUlrXdGLzBB/qV7DnfmIag789wvI9Vwq6LECClRBCCCGKoIplbFk3xpt+Xm4oCswPuljQJQESrIQQQghRRFlbmDG7TwNm9a6PhZnKsPzavYcFVpMEKyGEEEIUaQOaVWL1yOaG92VsLQusFglWQgghhCjy6lewN/xsb2NRYHVIsBJCCCGEMBEJVkIIIYQQJiLBSgghhBDCRCRYCSGEEEKYiAQrIYQQQggTkWAlhBBCCGEiEqyEEEIIIUxEgpUQQgghhIlIsBJCCCGEMBEJVkIIIYQQJiLBSgghhBDCRCRYCSGEEEKYiAQrIYQQQggTkWAlhBBCCGEiEqyEEEIIIUxEgpUQQgghhIlIsBJCCCGEMJEiF6xiYmIYOHAgGo0GBwcHRowYQUJCQrZ9kpOTCQgIoGzZstjZ2eHv709UVJTh87///psBAwbg5uaGjY0NderUYf78+U+sZ/fu3TRp0gQrKytq1KjBypUrTb17QgghhCjCilywGjhwIGfOnCEwMJAtW7awd+9eRo8enW2fyZMns3nzZtatW8eePXuIiIigd+/ehs9DQkJwdnZm9erVnDlzhrfffpvp06ezaNEiQ5urV6/SrVs32rdvT2hoKJMmTWLkyJFs3749z/ZVCCGEEEWLSlEUpaCLyK1z587h4eHB0aNH8fLyAmDbtm107dqVmzdv4urq+kSfuLg4nJycWLt2LX369AHg/Pnz1KlTh+DgYFq0aJHptgICAjh37hw7d+4EYOrUqfzxxx+cPn3a0KZ///7Exsaybdu2LGtOSUkhJSXF8D4+Ph43Nzfi4uLQaDRPfxCEEEII8YTE1HQ8ZugHO85+4IetpblJ1x8fH4+9vX2O399FasQqODgYBwcHQ6gC8PX1Ra1Wc/jw4Uz7hISEkJaWhq+vr2GZu7s7lSpVIjg4OMttxcXF4ejoaLTtf68DwM/PL9t1AMyaNQt7e3vDy83NLdv2QgghhCi6ilSwioyMxNnZ2WiZubk5jo6OREZGZtnH0tISBwcHo+UuLi5Z9jl48CA///yz0SnGyMhIXFxcnlhHfHw8SUlJWdY8ffp04uLiDK8bN25kt4tCCCGEKMIKRbCaNm0aKpUq29f58+fzpZbTp0/Ts2dP3nvvPTp16vTc67OyskKj0Ri9hBBCCFE8mfYE5DOaMmUKQ4cOzbZNtWrV0Gq1REdHGy1PT08nJiYGrVabaT+tVktqaiqxsbFGo1ZRUVFP9Dl79iwdO3Zk9OjRvPPOO0+s5993Ej5ah0ajwcbGJoc9FEIIIUResrU059qn3Qq6jMIRrJycnHBycsqxnbe3N7GxsYSEhODp6QnAzp070el0NG/ePNM+np6eWFhYEBQUhL+/PwBhYWGEh4fj7e1taHfmzBk6dOjAkCFD+PjjjzPd9tatW42WBQYGGq1DCCGEECVbkborEKBLly5ERUWxbNky0tLSGDZsGF5eXqxduxaAW7du0bFjR1atWkWzZs0AGDt2LFu3bmXlypVoNBomTJgA6K+lAv3pvw4dOuDn58fnn39u2JaZmZkh8F29epV69eoREBDA8OHD2blzJxMnTuSPP/7Az88v1/Xn9q4CIYQQQhQexfKuQIA1a9bg7u5Ox44d6dq1Kz4+PqxYscLweVpaGmFhYSQmJhqWzZ07l+7du+Pv70+bNm3QarWsX7/e8Pmvv/7KnTt3WL16NeXLlze8mjZtamhTtWpV/vjjDwIDA2nYsCFz5szh66+/fqpQJYQQQojirciNWBV1MmIlhBBCFD3FdsRKCCGEEKKwkmAlhBBCCGEiEqyEEEIIIUxEgpUQQgghhIlIsBJCCCGEMBEJVkIIIYQQJiLBSgghhBDCRCRYCSGEEEKYiAQrIYQQQggTkWAlhBBCCGEiEqyEEEIIIUzEvKALKGkePZoxPj6+gCsRQgghRG49+t7O6RHLEqzy2YMHDwBwc3Mr4EqEEEII8bQePHiAvb19lp+rlJyilzApnU5HREQEpUuXRqVS5em24uPjcXNz48aNG9k+iVs8HznO+UOOc/6RY50/5DjnD1MdZ0VRePDgAa6urqjVWV9JJSNW+UytVlOxYsV83aZGo5H/0+YDOc75Q45z/pFjnT/kOOcPUxzn7EaqHpGL14UQQgghTESClRBCCCGEiUiwKsasrKx47733sLKyKuhSijU5zvlDjnP+kWOdP+Q454/8Ps5y8boQQgghhInIiJUQQgghhIlIsBJCCCGEMBEJVkIIIYQQJiLBSgghhBDCRCRYFUG3bt3ilVdeoWzZstjY2FC/fn2OHTtm+FxRFGbMmEH58uWxsbHB19eXixcvGq0jJiaGgQMHotFocHBwYMSIESQkJOT3rhRq2R3ntLQ0pk6dSv369SlVqhSurq4MHjyYiIgIo3XIcc5ZTr/P/zZmzBhUKhXz5s0zWi7HOWe5Oc7nzp3jxRdfxN7enlKlStG0aVPCw8MNnycnJxMQEEDZsmWxs7PD39+fqKio/N6VQi+nY52QkMD48eOpWLEiNjY2eHh4sGzZMqN1yLHOXpUqVVCpVE+8AgICgNwdv/DwcLp164atrS3Ozs68+eabpKenP39xiihSYmJilMqVKytDhw5VDh8+rFy5ckXZvn27cunSJUObTz/9VLG3t1c2bNig/P3338qLL76oVK1aVUlKSjK06dy5s9KwYUPl0KFDyr59+5QaNWooAwYMKIhdKpRyOs6xsbGKr6+v8vPPPyvnz59XgoODlWbNmimenp5G65HjnL3c/D4/sn79eqVhw4aKq6urMnfuXKPP5DhnLzfH+dKlS4qjo6Py5ptvKsePH1cuXbqkbNy4UYmKijK0GTNmjOLm5qYEBQUpx44dU1q0aKG0bNmyIHap0MrNsR41apRSvXp1ZdeuXcrVq1eV5cuXK2ZmZsrGjRsNbeRYZy86Olq5ffu24RUYGKgAyq5duxRFyfn4paenK/Xq1VN8fX2VEydOKFu3blXKlSunTJ8+/blrk2BVxEydOlXx8fHJ8nOdTqdotVrl888/NyyLjY1VrKyslB9//FFRFEU5e/asAihHjx41tPnzzz8VlUql3Lp1K++KL0JyOs6ZOXLkiAIo169fVxRFjnNu5PY437x5U6lQoYJy+vRppXLlykbBSo5zznJznPv166e88sorWX4eGxurWFhYKOvWrTMsO3funAIowcHBJqu1qMvNsa5bt67ywQcfGC1r0qSJ8vbbbyuKIsf6Wbz22mtK9erVFZ1Ol6vjt3XrVkWtViuRkZGGNkuXLlU0Go2SkpLyXLXIqcAiZtOmTXh5efG///0PZ2dnGjduzFdffWX4/OrVq0RGRuLr62tYZm9vT/PmzQkODgYgODgYBwcHvLy8DG18fX1Rq9UcPnw4/3amEMvpOGcmLi4OlUqFg4MDIMc5N3JznHU6HYMGDeLNN9+kbt26T6xDjnPOcjrOOp2OP/74g1q1auHn54ezszPNmzdnw4YNhjYhISGkpaUZ/W1xd3enUqVKhr8tIne/0y1btmTTpk3cunULRVHYtWsXFy5coFOnToAc66eVmprK6tWrGT58OCqVKlfHLzg4mPr16+Pi4mJo4+fnR3x8PGfOnHmueiRYFTFXrlxh6dKl1KxZk+3btzN27FgmTpzI999/D0BkZCSA0S/Lo/ePPouMjMTZ2dnoc3NzcxwdHQ1tSrqcjvN/JScnM3XqVAYMGGB4yKcc55zl5jjPnj0bc3NzJk6cmOk65DjnLKfjHB0dTUJCAp9++imdO3fmr7/+olevXvTu3Zs9e/YA+uNsaWlp+IfDI//+2yJy9zu9cOFCPDw8qFixIpaWlnTu3JnFixfTpk0bQI7109qwYQOxsbEMHToUyN3xi4yMzPR78tFnz8P8uXqLfKfT6fDy8uKTTz4BoHHjxpw+fZply5YxZMiQAq6u+Hia45yWlkbfvn1RFIWlS5cWRLlFVk7HOSQkhPnz53P8+HFUKlUBV1t05XScdTodAD179mTy5MkANGrUiIMHD7Js2TLatm1bYLUXNbn527Fw4UIOHTrEpk2bqFy5Mnv37iUgIABXV1ejURaRO9988w1dunTB1dW1oEsBZMSqyClfvjweHh5Gy+rUqWO4c0er1QI8cfdDVFSU4TOtVkt0dLTR5+np6cTExBjalHQ5HedHHoWq69evExgYaBitAjnOuZHTcd63bx/R0dFUqlQJc3NzzM3NuX79OlOmTKFKlSqAHOfcyOk4lytXDnNz8xz/tqSmphIbG2vU5t9/W0TOxzopKYn/+7//48svv6RHjx40aNCA8ePH069fP7744gtAjvXTuH79Ojt27GDkyJGGZbk5flqtNtPvyUefPQ8JVkVMq1atCAsLM1p24cIFKleuDEDVqlXRarUEBQUZPo+Pj+fw4cN4e3sD4O3tTWxsLCEhIYY2O3fuRKfT0bx583zYi8Ivp+MMj0PVxYsX2bFjB2XLljVqL8c5Zzkd50GDBnHy5ElCQ0MNL1dXV9588022b98OyHHOjZyOs6WlJU2bNs22jaenJxYWFkZ/W8LCwggPDzf8bRE5H+u0tDTS0tJQq42/fs3MzAwjh3Ksc++7777D2dmZbt26GZbl5vh5e3tz6tQpo3+UPfrH8X+D8VN7rkvfRb47cuSIYm5urnz88cfKxYsXlTVr1ii2trbK6tWrDW0+/fRTxcHBQdm4caNy8uRJpWfPnplOt9C4cWPl8OHDyv79+5WaNWvK7en/ktNxTk1NVV588UWlYsWKSmhoqNFtv/++o0SOc/Zy8/v8X/+9K1BR5DjnJDfHef369YqFhYWyYsUK5eLFi8rChQsVMzMzZd++fYY2Y8aMUSpVqqTs3LlTOXbsmOLt7a14e3sXxC4VWrk51m3btlXq1q2r7Nq1S7ly5Yry3XffKdbW1sqSJUsMbeRY5ywjI0OpVKmSMnXq1Cc+y+n4PZpuoVOnTkpoaKiybds2xcnJSaZbKKk2b96s1KtXT7GyslLc3d2VFStWGH2u0+mUd999V3FxcVGsrKyUjh07KmFhYUZt7t27pwwYMECxs7NTNBqNMmzYMOXBgwf5uRuFXnbH+erVqwqQ6evRPCqKIsc5N3L6ff6vzIKVHOec5eY4f/PNN0qNGjUUa2trpWHDhsqGDRuMPk9KSlLGjRunlClTRrG1tVV69eql3L59O792ocjI6Vjfvn1bGTp0qOLq6qpYW1srtWvXVubMmaPodDpDGznWOdu+fbsCPPH9pii5O37Xrl1TunTpotjY2CjlypVTpkyZoqSlpT13XSpFUZTnG/MSQgghhBAg11gJIYQQQpiMBCshhBBCCBORYCWEEEIIYSISrIQQQgghTESClRBCCCGEiUiwEkIIIYQwEQlWQgghhBAmIsFKCCGEEMJEJFgJIYQQQpiIBCshRIFQqVTZvmbOnFnQJZpclSpVmDdvXkGX8YRLly4xbNgwKlasiJWVFVWrVmXAgAEcO3asoEsTosgxL+gChBAl0+3btw0///zzz8yYMYOwsDDDMjs7u4Io66kpikJGRgbm5vn35zQ1NRVLS0uTrOvYsWN07NiRevXqsXz5ctzd3Xnw4AEbN25kypQp7NmzxyTbEaKkkBErIUSB0Gq1hpe9vT0qlcpo2U8//USdOnWwtrbG3d2dJUuWGPpeu3YNlUrFL7/8QuvWrbGxsaFp06ZcuHCBo0eP4uXlhZ2dHV26dOHOnTuGfkOHDuWll17i/fffx8nJCY1Gw5gxY0hNTTW00el0zJo1i6pVq2JjY0PDhg359ddfDZ/v3r0blUrFn3/+iaenJ1ZWVuzfv5/Lly/Ts2dPXFxcsLOzo2nTpuzYscPQr127dly/fp3JkycbRuUAZs6cSaNGjYyOzbx586hSpcoTdX/88ce4urpSu3ZtAG7cuEHfvn1xcHDA0dGRnj17cu3atVz/N1AUhaFDh1KzZk327dtHt27dqF69Oo0aNeK9995j48aNuV6XEEJPRqyEEIXOmjVrmDFjBosWLaJx48acOHGCUaNGUapUKYYMGWJo99577zFv3jwqVarE8OHDefnllyldujTz58/H1taWvn37MmPGDJYuXWroExQUhLW1Nbt37+batWsMGzaMsmXL8vHHHwMwa9YsVq9ezbJly6hZsyZ79+7llVdewcnJibZt2xrWM23aNL744guqVatGmTJluHHjBl27duXjjz/GysqKVatW0aNHD8LCwqhUqRLr16+nYcOGjB49mlGjRj31MQkKCkKj0RAYGAhAWloafn5+eHt7s2/fPszNzfnoo4/o3LkzJ0+ezNWIVmhoKGfOnGHt2rWo1U/+O9vBweGp6xSixFOEEKKAfffdd4q9vb3hffXq1ZW1a9catfnwww8Vb29vRVEU5erVqwqgfP3114bPf/zxRwVQgoKCDMtmzZql1K5d2/B+yJAhiqOjo/Lw4UPDsqVLlyp2dnZKRkaGkpycrNja2ioHDx402vaIESOUAQMGKIqiKLt27VIAZcOGDTnuV926dZWFCxca3leuXFmZO3euUZv33ntPadiwodGyuXPnKpUrVzaq28XFRUlJSTEs++GHH5TatWsrOp3OsCwlJUWxsbFRtm/fnmNtiqIoP//8swIox48fz1V7IUTOZMRKCFGoPHz4kMuXLzNixAijkZ309HTs7e2N2jZo0MDws4uLCwD169c3WhYdHW3Up2HDhtja2hree3t7k5CQwI0bN0hISCAxMZEXXnjBqE9qaiqN/7+duwlJrIvDAP5klFgZNGRBURRpcAOpsSgIqkWL2qQUgZAbwdpFELiNvkghykVhmwJz0SKCaKMt2kQRRERokYSkhRQFtXAhKH05i6GGm+/bxzu+MzLz/Faeew/nnOvq4dz/PV+/iq7V1dWJ2pFIBCMjI3C5XLi6usLDwwOi0ShCodBnHv9fqdVq0S6U1+vF6ekp5HK5qF8sFkMgEPjQmPF4PClrI6IfGKyIKKVEIhEAwPz8PBoaGkT30tPTRe2MjIyX3881S6+vPT09fXpul8uF4uJi0T2pVCpqZ2dni9pmsxkbGxuYmpqCUqmETCZDd3e3qH7rn0gkkoSAc39/n9Dv9XyRSAS1tbVYWlpK6KtQKN6c81llZSUA4OTkJCE4EtF/w2BFRCmlsLAQRUVFCAaDMBgMSR/f6/UiGo1CJpMBAHZ3d5GTk4OSkhJ8+fIFUqkUoVBIVE/1ETs7OzAajejs7ATwPfi8LiTPzMzE4+Oj6JpCocD19TXi8fhLOPR4PO/Op9FosLy8jIKCAuTm5n5qrc9qampQVVWF6elp6PX6hDqrcDjMOiuiT+JXgUSUckZHR2G1WjEzMwO/34+joyM4HA7YbLafHvvu7g4mkwk+nw9utxvDw8Po7++HRCKBXC6H2WzG4OAgnE4nAoEADg4OMDs7C6fT+ea4KpUKq6ur8Hg88Hq96OnpSdgtKysrw9bWFi4vL3F7ewvg+9eCNzc3mJycRCAQgN1ux/r6+rvPYTAYkJ+fD51Oh+3tbZydnWFzcxMDAwO4uLj40H+RlpYGh8MBv9+PpqYmuN1uBINBHB4eYmJiAjqd7kPjENEPDFZElHJ6e3uxsLAAh8MBtVqNlpYWLC4uory8/KfHbm1thUqlQnNzM/R6PbRaregw0vHxcQwNDcFqtUIQBLS3t8Plcr07t81mQ15eHhobG9HR0YG2tjZoNBpRn7GxMZyfn6OiouLldZ0gCJibm4Pdbkd1dTX29vZgNpvffY6srCxsbW2htLQUXV1dEAQBJpMJsVjsZQfr+WiIt45gqK+vx/7+PpRKJfr6+iAIArRaLY6Pj1PyMFOiVJcWZ/UiEf0ljEYjwuEw1tbWfvdSfgmHwwGLxQKfzyeqPSOi/w93rIiI/lButxsWi4WhiugXYvE6EdEfamVl5Xcvgeivw1eBREREREnCV4FEREREScJgRURERJQkDFZEREREScJgRURERJQkDFZEREREScJgRURERJQkDFZEREREScJgRURERJQk3wBClJjobJTDQAAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 640x480 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "# Calculate \n",
    "import numpy as np\n",
    "from scipy.optimize import curve_fit\n",
    "import matplotlib.pyplot as plt\n",
    "\n",
    "# Define linear function\n",
    "def func (x,a,b):\n",
    "    return a*x +b\n",
    "\n",
    "x = np.array(temperatures)\n",
    "# Convert keff into reactivity\n",
    "y = np.array([(x[0]-1)/x[0] for x in keff])\n",
    "yerr = np.array([x[1] for x in keff])\n",
    "\n",
    "#Fit parameters of the function to the data\n",
    "popt, pcov = curve_fit(func, x, y)\n",
    "\n",
    "# Calculate average isothermal coefficient in pcm/C \n",
    "alpha_iso = popt[0]*1e5\n",
    "\n",
    "plt.figure()\n",
    "plt.errorbar(x, y, yerr=yerr, label='Simulation data')\n",
    "plt.plot(x, func(x, *popt), '--', label='{:.3f} [pcm/C]'.format(alpha_iso))\n",
    "plt.legend()\n",
    "plt.xlabel('Temperature, C')\n",
    "plt.ylabel(r'$\\rho\\,\\pm\\sigma$')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "da6d7a55",
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3 (ipykernel)",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.10.7"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}
